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Classical  thermodynamic  concepts  were  used  to  derive  general  rela- 
tionships to  describe  the  equilibrium  disposition  of  charged  species  in 
aqueous  contact  with  a  cation  exchanger.   It  vjas  shown  through  these 
relationships   that  anions  would  be  excluded  by  the  exchanger  and  solu- 
tion phase  cations  of  the  same  type  as  those  initially  saturating  the 
exchanger  influenced  the  adsorption  of  counterions .   The  latter  result 
was  verified  vrith  exchange  adsorption  isotherms  for  Na"*"  and  Li   in  solu- 
tions of  0.05  M,  0.02  M,  and  0.005  M  CaCNO^)^  and  a  Ca-saturated  ex- 
changer.  These  results  illustrated  that  Na'^  and  Li   adsorption  de- 
creased v;ith  increasing  concentrations  of  Ca^^  in  the  equilibrium  solu- 
tion.  Similar  results  were  observed  for  the  adsorption  of   Ca   in 
solutions  of  0.07  5  M  and  0.05  M  C.^^(^0    )^.      The  adsorption  isotherms  were 
linear  over  the  range  of  concentrations  studied  for  a  given  Ca(NO^) 
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concentration.   The  thermodynamic  relations  in  conjunction  with  the 
Dcbve-Huckel  theory  were  used  to  predict  the  slopes  for  these  isotherms. 
The  computed  values  agreed  reasonably  well  with  their  experimental 
counterparts.   These  results  demonstrated  inductively,  that  electro- 
static ion-ion  interactions  as  conceived  and  quantified  by  the  Debye- 
Huckel  theory,  constituted  the  physical  basis  fox  the  observed  equi- 
librium  behaviour. 

The  consequences  of  the  foregoing  results  on  the  transport  be- 
haviour of  Na+,  Li+,  '*5ca2+  and  CI"  were  investigated  using  miscible 
displacenient  experiments  in  laboratory  columns  packed  with  the  exchanger 
material.   These  studies  were  conducted  using  pulse  inputs  of  tracer 
solutions  containing  various  concentrations  of  Ca(N03)2.   Tritiated 
water  (HTO)  was  introduced  with  all  tracer  pulses  to  evaluate  hydro- 
dynamic  dispersion.   An  analytical  solution  to  the  convective-dispersion 
mass  transport  equation  for  a  reactive  solute  (linear  adsorption  iso- 
thern)  was  used  to  describe  the  experimental  data. 

A3  expected,  the  elution  curves  for  CI'  were  displaced  to  the  left 
of  those  for  HTO  illustrating  exclusion  of  CI"  by  the  exchanger.   The 
elution  curves  for  Na+,  Li+  and  ^Sca^*-  were  displaced  to  the  right  of 
those  for  HTO,  the  shift  increasing  with  decreasing  concentration  of 
C;a2+  in  the  tracer  solution.   Anomalous  and  unusual  patterns  were  ob- 
served in  some  of  the  elution  curves.   The  elution  curves  for  Na+  in 
0.005  M  CaCNO^)^  and  for  ^5ca2+  in  0.05  M  Ca(N03)2  showed  deviations 
that  simulated  the  effect  of  kinetic  mass  transfer  processes.   Inflec- 
tions were  observed  on  the  desorption  side  of  the  pulse  elution  data 
for  Li+  in  0.02  U   and  0.005  M  CaCNO^)^  and  for  ^5ca2+  m  0.075  M 
Ca(XO  )  .   These  elution  curves  were  not  described  by  the  analytical 


solution.   Further  experimentation  demonstrated  that  these  data  could 
be  explained  on  the  basis  of  differences  between  the  parameters 
characterizing  exchange  adsorption  and  desorption.   A  quantitative 
treatment,  based  on  the  thermodynamic  relations  coupled  with  the 
Debye-Huckel  theory,  showed  that  these  differences  were  a  direct  con- 
sequence of  the  metathetical  nature  of  the  sorption  process. 

Some  cursory  observations  were  made  regarding  the  transport  be- 
iiaN/ior  of  Na  ,  1!T0  and  Cl~  under  steady-state  unsaturated  water  flow 
conditions.   The  exclusion  of  Cl~  was  again  evident.   Discrepancies 
between  the  experimental  and  analytical  curves  for  Na  indicated  that 
portions  of  the  water-unsaturated  exchanger  phase  became  either  inacces- 
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INTRODUCTION 

The  cla.sjuc  work  of  Thonpson  and  IJay,  repoi-tecl  during  1850-1852 
en  exchange  reactlono  involving  soil  materials  treated  with  inorganic 
salt  solutions,  marked  the  beginning  of  syster.iatic  investigations  of 
ion  ex.--hr.nge.   Their  experinents  were  conducted  decades  before  Arrenhius 
proposed  the  theory  of  electrolytic  dissociation  and  r.'.any  years  before 
the  law  of  mass  action  was  enunciated.   There  was,  therefore,  no  exist- 
ing theoretical  framework  within  which  the  tv.'o  English  agricultural 
chemists  could  interpret  their  observations.   The  principle  that  sub- 
stances did  not  react  except  in  a  dissolved  state  v/as  generally  accepted 
by  chemists  of  their  era.   Ifnether  the  reactions  reported  were  purely 
chemical  or  physical  became  a  debated  question.   This  dichotomy  in 
views  has  still  not  been  fully  resolved. 

For  a  period  after  its  discovery,  ion  exchange  remained  of 
academic  interest  until  about  1907  when  it  was  used  to  soften  water. 
Instability  under  non-neutral  pH  conditions  and  the  lov;  capacity  of  the 
available  natural  and  artificial  siliceous  ion  exchangers  provided  the 
impetus  for  the  discovery  and  synthesis  of  organic  ion  exchangers.   These 
synthetic  materials  are  now  utilized  in  most  industrial  and  laboratory 
applications  involving  ion  exchange.   Extensive  research  concurrent 
with  their  increasingly  widespread  use  has  resulted  in  the  evolution 
of  more  comprehensive  and  refined  ion  exchange  theories. 

The  miscible  displacement  technique  has  in  the  past  few  decades 
become  an  increasingly  valuable  experiiaental  tool  for  studying 
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nh">-3ico-c)ic'iiical  processes  ai^sociated  with  solute  transport  in  porous 
rnoiJia.   Theoretical  concepts  on  physico-chemical  interactions  of  partic- 
ular solutes  have  been  successfully  incorporated  into  classical  mass 
tra-i..,port  models.   Comparisons  between  experimental  and  theoretical 
results  on  solute  transport  have  served  to  verify  concepts  and  to 
evaluate  the  Influence  of  physical  factors.   The  diversity  of  possible 
physico-chemical  interactions  of  inorganic  ions  v.'ith  various  soil 
components  makes  natural  soil  materials  a  complex  medium  in  which  to 
study  ion  exchange.   This  difficulty  can  be  circumvented  somewhat  by 
the  use  of  synthetic  exchangers. 

The  primary  objective  of  this  study  was  to  evaluate  some  theore- 
tical concepts  on  ion  exchange  equilibria  and  kinetics  using  the  mis- 
cible  displacement  technique  and  a  porous  medium  prepared  from  a  syn- 
thetic organic  ion  exchanger. 


CIRPTER  1 

A  REVIEW  OF  PERTINENT  COKCEPTS  ON  SOLUTE  TPJVNSPORT  DURING  KISCIBLE 

DISPLACEMENT  IN  DISCRETE  POROUS  MEDIA  AND  ON  INORGANIC  ION  EXCHAIn'GE 

EQUILIBRIA  AND  KINETICS 

Miscihle  Displacement  Processes  in  Discrfcte 
Porous  Media 

Miscible  displacement  is  the  term  used  to  describe  a  process 
whereby  one  fluid  is  displaced  by  another,  both  fluids  being  miscible 
in  each  other.   This  process  occurs  in  a  soil  when  a  solution  is  dis- 
placed downwards  by  incoming  rain  or  irrigation  water  or  v/hen  sea  water 
displaces  fresh  water  during  drawdoum  in  coastal  aquifers.   In  the 
petroleum  industry,  fluids  miscible  with  crude  oil  are  used  to  increase 
the  efficiency  with  which  the  crude  oil  is  displaced  from  oil  bearing 
strata  during  secondary  operations.   Separation  and  recovery  processes 
in  industrial  packed  towers  or  in  chromatographic  colun-ius  are  other 
examples  of  miscible  displacement  processes  in  porous  media. 

The  transport  behavior  of  materials  dissolved  in  the  displacing 
or  displaced  fluids  is  influenced  by  the  physical  properties  of  the 
medium  (particle  size,  shape,  and  m? vner  of  packing),  rheological  pro- 
perties of  the  fluids,  and  both  the  equilibria  and  rates  of  reactions 
between  the  solute  and  solid  matrix  of  the  porous  medium.   Miscible 
displacement  experiments  where  there  is  no  interaction  between  the 
medium  and  the  solute  can  be  used  to  evaluate  the  effect  of  the  first 
two  factors.   The  combined  effect  of  the  first  two  factors  has  been 
grouped  under  the  term  "dispersion"  and  a  considerable  body  of  literature 
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has  gro\\Ti  out  of  these  studies.   Excellent  reviev;s  on  dispersion  have 
been  presented  by  Fried  and  Combarnous  (25)  and  Bear  (6). 

Dispersion  of  Solutes  During  Miscible 
Displacement  in  Discrete  Porous  Media 

Dispersion  is  the  result  of  a  physical  mixing  between  displacing 
and  displaced  fluids.   Its  effect  can  be  studied  by  measuring  the  rela- 
tive concentration  gradients  between  the  displacing  and  displaced 
fluids.   Gradients  are  usually  determined  experimentally  by  withdrawing 
samples  at  suitable  time  intervals  at  a  fixed  location  in  the  flov;  region 
under  consideration.   This  technique  and  the  relative  concentration 
versus  time  plots,  termed  breakthrough  curves,  have  been  described  by 
Nielsen  and  Biggar  (53). 

Dispersion  in  many  practical  situations  involving  solute  transport, 
especially  in  chromatographic  and  industrial  separations,  has  provided 
much  of  the  initial  motivation  to  quantify  and  describe  the  process. 
Sevetal  theoretical  approaches  have  been  taken  in  this  connection. 
Attempts  have  been  made  to  predict  macroscopic  dispersion  effects  based 
on  various  geometrical  models  of  the  microstructure  of  the  channels  in 
the  porous  matrix  (23,33,64).   These,  however,  require  complex  mathe- 
matical analyses,  are  rather  academic  in  nature,  and  have  not  proven  to 
be  useful.   Limitations  of  this  approach  have  been  discussed  by  Scheidegger 
(66).   Analysis  of  solute  transport  based  on  a  physical  description  of 
the  flow  region  into  theoretical  plates  has  proved  quite  useful  in  eval- 
uating dispersion  in  chromatographic  separations  (28,68). 

An  approach  based  on  differential  mass-balance  equations  for  the 
solute  in  the  fluid  flow  field  has  been  successful  and  has  gained  wide 


accf-pLonce.   Taylor  (67)  used  this  approacli  to  study  solute  dispersion 
during  laminar  flow  in  straight  capillary  tubes.   Analysis  of  the 
differential  mass-balance  equation  with  appropriate  initial  and  boundary 
conditions  yielded  expressions  which  satisfactorily  explained  his  expori- 
m.-ntal  results.   These  results  shov;ed  that  dispersion  v/as  dependent  on 
the  fluid  velocity  and  the  distance  travelled  by  the  displacing  front. 
During  laminar  flov/  of  Newtonian  fluids  in  straight  tubes, a  typical 
parabolic  velocity  distribution  is  produced.   This  would  not  be  the  case 
for  a  discrete  porous  medium  with  geometrically  complex,  tortuous  channels 
bet<-.eon  the  pai:ticles.   For  a  porous  medium  it  is  assumed  that  the  micro- 
scopic velocities  between  particles  would  fluctuate  continuously  in  a 
random  manner  about  some  mean  value.   This  viewpoint  led  to  a  concept 
that  \v7as  developed  almost  simultaneously  by  Scheidegger  (65)  and 
Danckwerts  (16).   Although  they  used  somewhat  different  approaches  in 
their  arguments,  they  agreed  that  dispersion  could  be  regarded  as  a 
quasi-dif fusion  process.   With  this  concept,  the  classical  differential 
eciuation  for  mass  transport  by  simultaneous  coiivection  and  diffusion 
could  be  applied  if  an  appropriate  change  was  made  in  the  physical 
meaning  attached  to  the  diffusion  coefficient.   The  validity  of  tliis 
approach  was  demonstrated  immediately  by  other  investigators  (4,19,52). 

The  Convective-Dispersion  Equation  Including 

Assumptions,  Auxiliary  Conditions,  Analysis , 

and  Limitations 

Mass  balance  considerations  for  transport  by  simultaneous  diffu- 
sion and  convection  in  a  homogenous,  isotropic  porous  medium  yield 
the  classical  convective-dif fusion  equation.   In  one-dimension  this 
equation  (39)  is 


vC}  (1) 

3t    3x    8x     - 

where,  C(x,t)  is  concentration  of  solute  in  the  fluid  of  the  flov/  re- 
gion (M/l3),  t  is  time  (T) ,  x  is  distance  (L) ,  D  is  the  diffusion 
coefficient  (L^/T),  and  v  is  vector  velocity  of  the  fluid  at  any  point 
in  the  flow  region  (L/T) . 

Because  of  the  pore  size  distribution,  complex  fluid  velocity 
fields  arise  within  the  channels  of  the  porous  medium.   It  is  assumed, 
however,  that  these  velocity  vectors  are  statistically  distributed 
about  some  fixed  mean  velocity  vector  v* .   The  probability  distribu- 
tion of  the  deviations  about  this  average  is  then  considered  by  the 
introduction  of  a  quasi-dif fusion  coefficient,  D* ,  which  characterizes 
the  unidirectional  transport  of  the  solute.   This  quasi-dif fusion  co- 
efficient has  been  termed  the  "dispersion  coefficient"  by  Schcidegger 
(65).   If  it  is  also  assumed  that  the  dispersion  coefficient  is  inde- 
pendent of  concentration,  then  equation  (1)  reduces  to 

IC  ^  D*^  -  v*lC  (2) 

As  shown  by  Scheidegger  (65)  and  in  a  somewhat  more  heuristic  manner 
by  Rifai  et  al.  (59),  this  equation  can  also  bo  obtained  by  arguments 
based  entirely  on  probability  calculus.   This  equation  has  been  appro- 
priately called  the  "convective-dispersion"  equation. 

It  is  at  once  apparent  from  its  description  that  v''  is  the  ex- 
perimentally measurable  average  pore-water  velocity  of  the  fluid  in 
the  porous  medium.   This  is  given  by  the  relationship 


^*  =67^  (3) 

3 
where,  Q  is  volume  of  fluid  introduced  (L  ) ,  A  is  the  cross  sectional 

2 
area  of  the  flow  region  (L  ) ,  6   is  the  fractional  fluid  filled  poro- 

3   3 
sity  of  tne  medium  (L  /L  ),  and  t  is  time  (T) . 

It  is  also  clear  that  D  is  related  to  D"-"'  in  the  following  manner: 

—  '    V"  ->  0  (4) 

For  a  fixed  X  =  L,  the  solution  of  the  convective-dispersion 
equation  would  be  C  =  C(t,  D*,  L,  v*)  where  D*,  L,  and  v*  are  para- 
meters „   It  is  advantageous  to  formulate  the  differential  equation 
and  auxiliary  conditions  in  dimensionless  form.   The  problem  then 
becomes  independent  of  units  and  the  parameters  are  reduced  to  dimen- 
sionless groupings.   In  addition,  interactions  among  the  parameters 

become  more  obvious.   Tlie  dimensionless  variables  are   C*  =  -^, 

Co 

X  =  ^,  and  T  =  ^^^  which,  when  substituted  into  equation  (2)  give 
L  L 

dT  v*L     J^  8X   .  (5) 

Tlic  reciprocal  of  the  dimensionless  grouping  —y:     has  been  termed 
the  Peclet  number,  usually  denoted  by  P,  and  appears  as  a  single 
parameter  in  solutions  of  the  transformed  equation.   The  dimension- 
less time  T  corresponds  physically  to  the  number  of  fluid-filled 
void  volumes  introduced  into  the  flow  region. 

The  convective-dispersion  equation  in  either  form  [equation 
(2)  or  (5)]  has  been  analyzed  for  various  sets  of  boundary  and  ini- 
tial conditions  appropriate  to  different  physical  situations.   These 
have  been  summarized  by  Nir  and  Gershon  (54).   Tlie  analyses  are  based 
on  the  theory  of  equations  in  partial  derivatives,  a  very  difficult 


and  still  Incomplete  branch  oi'  higher  mathematics.   Analytical  solu- 
tions can  be  obtained  with  comparative  ease  for  simple  sets  of  boundary 
and  initial  conditions;  however,  for  complex  conditions  or  for  finite 
domains  of  definition,  one  is  forced  to  resort  to  numerical  integration 
methods. 

Analytic  solutions  to  equation  (2)  have  been  obtained  for  two 
simple  sets  of  conditions.   These  can  be  written  as 

I.   C(x  >  0,  t  =  0)  =  0     ,     C(x  =  0,  t  >  0)  =  Co, 

C(x  -  ~,  t  >  0)  =  0  ^^^ 

II.   C(x  >  0,  t  =  0)  =  0    ,     C(x  =  0,  0  <  t  ^  tj)  =  Co 

C(x  =  0,  t  >  t^)  =  0    ,     C(x  =  <-,  t  >  0)  =  0      ^■^^ 

Physically,  tlie  first  set  of  conditions  corresponds  to  a  solute 
of  concentration  Co  continuously  displacing  the  solvent  in  a  semi- 
infinite  medium.   In  the  second  case,  the  same  displacement  is  allowed 
to  proceed  for  a  time  t   and  for  all  times  thereafter  the  displacing 
solution  is  replaced  by  pure  solvent.   These  two  situations  have  been 
termed  miscible  displacement  with  a  "step  input"  and  "pulse  input," 
respectively. 

The  solution  to  equation  (2)  for  a  "step  input"  (59)  is 

given  by 

C(x,t)      If   ,  r>^-vt-,   ,      ,vx.    ,   x+vt-.  ■■         .  , 
-to  =-{erfc[^;^]   +  exp(-)  erfcl-,-^-]}         (8) 

where  the  superscripts  have  been  omitted  from  D"-''  and  v"'-'. 

Except  near  the  inlet  (x  =  0)  or  for  very  small  times,  unless 
D  is  large,  the  second  term  in  equation  (8)  is  small  and  can  be 
neglected.   Applying  this  condition  and  the  identity  erfc(x)  = 
1-  erf(x),  equation  (8)  reduces  to 


The   cumulative   distribution   function   r(;:)    for   a   standardized 
Gaussian   random  variable   of   mean   x  and   standard   deviation   a   is   jp'ven 
by 

P(x)      =     ^{1  +  erf[^2]>    °         .  (10> 

From  the  identity  erf(-x)  =  -erf(x),  it  is  apparent  that  equa- 
tion (9)  represents,  for  a  given  value  of  x  =  L,  a  Gaussian  cumula- 
tive distribution  function  of  mean  L  and  standard  deviation  /2Dt. 

The  solution  for  the  "pulse  input"  is  obtained  by  shifting  the 
solution  for  the  "step  input"  by  the  time  width  of  the  pulse  and  sub- 
tracting this  from  the  original  solution.   Applying  this  ratlier  in- 
tuitive procedure  to  equation  (9)  gives 

where  U(t-t  )  is  the  Heaviside  unit  function. 

The  previous  equations  have  been  used  to  analyze  experimental 
results  from  finite  laboratory  columns.   As  shown  by  Nir  and  Gershon 
(54)  only  a  small  error  is  involved  by  assuming  a  finite  rather  than 
a  semi-infinite  column. 

The  boundary  conditions  appropriate  for  a  defined  finite  domain 
arc  obtained  by  imposing  the  conservation  of  mass  on  the  fluxes  at  the 
boundaries.   This  results  in  the  following  boundary  conditions  for  a 
step  input  and  a  space  domain  of  definition  0  <^  x  <^  L 

x  =  0;t>0;    vCo  =  vC|^^Q-I>f^:^^0  (12a) 


.3C 


x=L;t>0;    vCe  =  vC|  ^  ^  ^  -  1>^!  _  ^  (12b) 

whore  Ce  is  the  exit  concentration  at  a  given  time. 


Ho\;ever,  as  argued  ])y  Danckv;erts  (16),  if  at  x  =  L,  "^1^^=^  "^  ^ 
thtn  C|  _   <  cl  ^  ^    A  sinilar  argunent  holds  if  at  x  =  L, 

^—l  >   0.   Hence  -;--   ,  must  be  equal  to  zero  otherwise  a  maximum 

■c)x|x=h  9x1  x=L 

or  a  niinimuTP,  exists  in  the  interior  of  the  column.   The  appropriate 
condition  at  x  =  L  must  therefore  be  -~\    ^,  =  0. 

An  approximate  solution  to  the  convective-dispersion  equation 
ill  dimensionless  form  for  the  above  boundary  conditions  has  been  pre- 
sented by  Brenner  (l-'t).   In  dimensionless  form  these  conditions  become 

X  =0;   T   >  0;    C'^  -  ~~     =  1  (13a) 

X  =  1;    T   >  0;   |f^  =  0  (13b) 

Wehner  and  Uilhelm  (75)  and  Lindstrom  et  al.  (46)  have  shomi  that, 
for  large  values  of  P,  the  condition  at  x  =  0  reduces  to  C"  =  1  for 
t  >  0. 

In  order  to  make  the  condition  at  x  =  0  homogeneous,  Brenner 

defines  C-  =  77-"^  where  C.  is  the  initial  concentration  in  the 

'i 
fluid  at  t  =  0.   The  initial  concentration  was  taken  to  be  P'.ero 

for  his  solution.   Brenner  tabulated  numerical  values  of  his  solu- 
tion for  a  wide  range  of  values  of  P/4.   It  was  found  that,  v;ith 
increasing  values  of  P,  his  solution  asyraptotically  apj^roached 
the  simplified  solution  given  by  equation  (9).   He  also  discusses  the 
limiting  behavior  of  his  solution  for  P  =  0  and  P  =  ".   These  corre- 
spond physical ly  to  complete  instantaneous  mixing  and  to  no  dispersion, 
respectively.   His  solution  has  been  used  frequently  for  analysis  of 
experimental  data  from  misciblo  displacement  experiments  (61). 

Numerical  procedures  for  integration  of  the  convective-dispersion 
equation  are  all  based  on  methods  whereby  continuous  systems  are 
matlienatically  reduced  to  equivalent  discrete  systems.   Suitable  algo- 


rithms  can  then  be  developed  and  solved  itcrativcly  using  high  speed 
computers.   Serious  and  difficult  mathematical  questions  of  adequacy, 
accuracy,  convergence  and  stability  arise  for  numerical  solutions  and 
adequate  answers  are  often  not  available.   Procedures  based  on  finite 
differences  and  finite  elements  are  discussed  by  Ames  (1)  and  Finder 
and  Gray  (55) . 

From  the  assumptions  used  in  its  development,  the  convective- 
dispersion  equation  represents  an  idealized  conception  of  the  miscible 
displacement  process.   Its  shortcomings  were  demonstrated  early  by 
its  inability  to  describe  non-reactive  solute  behavior  in  unsaturated 
porous  media  (53).   Also,  it  fails  to  predict  the  effects  of  density 
differences  between  displacing  and  displaced  fluids  (62).   Its  major 
weakness  probably  lies  in  the  assumption  of  a  fluid  continuum.   The 
presence  of  dead-end  pores  and  regions  where  the  fluid  is  hydro- 
dynamically  immobile  would  clearly  lead  to  unpredicted  results.   R.e-- 
finements  made  to  handle  such  cases  (15,69)  have  resulted  in  better 
agreement  between  predicted  and  experimental  results. 

Extension  of  the  Continuum  Mass 


Balance  Approach  for  Single  Interacting 
Solutes 


The  convective-dispersion  equation  can  be  extended  to  describe 
the  transport  of  a  single  interacting  solute  through  a  porous  medium. 
This  requires  additional  terms  to  cover  the  time  rate  of  change  of 
solute  concentration  in  the  fluid  as  a  result  of  mass  transfer  be- 
tween the  fluid  and  solid  matrix  of  the  medium.   For  p.   general  case, 
these  terms  constitute  the  difference  between  the  instantaneous  adsorp- 
tion and  desorption  rates.   Equivalently  and  more  conveniently,  this 


can  be  expre5;rfed  in  teriiis  of  the  instantaneous  net  time  rate  of 
accumulation  of  the  solute  on  the  particles  of  the  raediuni.   The 
concentration  of  the  solute  in  the  matrix  expressed  as  mass  of 
solute/unit  mass  of  matrix  can  be  denoted  by  S(x,t).   The  instan- 
taneous net  time  rate  of  change  of  concentration  in  the  fluid  due 
to  mass  transfer  is  then  -r—  ——   in  which  p  is  the  dry  bulk  density 
of  the  matrix  and  6   is  the  fractional  water-filled  porosity.   Con- 
sidering only  one  mass  transfer  process,  the  resulting  differential 
equation  is 

i£  =  n-^^c   _   ...ic  _  P__  is 

9t       -3^      ^'  3x     e^  9t  (14) 

The  functional  form  of  S(x,t)  has  profound  consequences  on  the  analysis 
of  the  resulting  differential  equation.   In  general,  the  mass  transfer 
processes  may  involve  purely  diffusion  kinetics,  irreversible  and 
reversible  chemical  kinetics  or  both.   Hence,  it  would  be  expected 
that  in  the  general  case 

~     =      f[C,  S,  T,  (A,,  A,,. ..A  ),..r      .       (A,,  A^...A,,)  ,   .   J   ,^,. 
dt  '  ''  ^  1'   2     n  diffusxon,  1    2     ^1  chemical    (15) 

where  1    is  temperature,  and  A  ...A  are  parameters  characteris;ing 

diffusion  and  chemical  kinetics  (57).   The  case  of  mass  transfer 

involving  only  irreversible  chemical  kinetics  has  been  analyzed  by 

Anundson  (2). 

Processes  involving  only  reversible  chemical  kinetics  can  be 

represented  in  the  general  case  as 

3S 


3t 


A^f^(S,  C,  T  I   -   A2f2(^^  C,  t)  (16) 


where  A,  and  A   are  constants  characterizing  the  kinetics  of  sorp- 
tion and  dcsoi  ption,  resjiectively . 
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An  important  simplification  of  chemical  kinetics  is  based  on  tVie 
assumption  that,  at  a  constant  temperature,  the  processes  occur  fast 
enough  to  insure  that  equilibrium  is  instantaneous.   The  functional 
form  of  S(x,t)  is  then  given  by  an  equation  which  describes  the  equi- 
librium isotherm.   This  equation  can  be  experimentally  determined 
or  deduced  from  kinetic  equations,  if  these  are  knovni^  by  setting 

— —  =0.   In  general,  these  isotherms  can  be  represented  as 
a  t 


=   f(C,  A^...A^)  (17) 


where  X    ...X     are  parameters  characterizing  the  equilibrium  distri- 
bution of  the  solute  between  the  fluid  and  solid  matrix  of  the  medium. 

at     dt        dt 


With  this  assumption,  introducing  -r—  =  TT  °  TT  in  the  differen- 


tial equation  gives 

_9C  ^  D2_  3^C     V"  9C 

at     R  J^     ~     R   9x  (18) 

where,  R  =   1  +  ^  -  .ii 

From  this  equation,  some  useful  deductions  can  be  made  regarding 
the  transport  behavior  of  the  solute  for  certain  generalized  func- 
tional forms  for  the  isotherm. 

In  the  case  v/here  the  isotherm  is  a  single-valued  concave  dovm- 

j  ^   .   osl     as  I     .  .   v'>|     "v-t   ^   ^  ^  „ 

ward  function  -; —     >  -r—     and  thus  -r—     <  •=—    tor  C,  <   L„. 
ctc|cj^     9c|c2  ^    l^-,  '^    >^2 

Similarly  the  reverse  is  true  for  C^  <  C  v;hen  an  adsorption  isotherm  is 
represented  by  a  concave  upward  function.   For  a  linear  isotherm, 
— -  is  independent  of  C.   Typically,  for  a  pulse  input  of  a  non-reactive 
solute,  dispersion  causes  the  elution  curves  to  have  a  gradual  increase 
in  concentration  from  zero  to  a  maximum  on  the  front  and  falling  again 
to  zero  on  the  back  side  of  the  elution  data.   For  a  concave  downward 
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isotherm,  the  concentration  effect  on  the  velocity  would  result  in 
a  nullification  of  dispersion  on  the  front  and  an  enhancement  on  the 
back  side  of  the  elution  curve.  This  effect  would  not  appear  in  the 
case  of  a  linear  isotherm  and  the  reverse  would  be  true  for  a  concave 
upward  isotherm.  For  a  linear  isotherm,  the  front  and  rear  of  the 
elution  curve  would  match  exactly,  but  will  not  match  for  non-linear 
isotherms.  If  the  assumption  of  instantaneous  equilibrium  was  valid 
for  a  particular  situation,  the  shape  of  the  elution  curve  can  pro- 
vide a  valuable  indication  of  the  type  of  adsorption  isotherm  in- 
volved in  the  displacement  process.  These  arguments  have  been  con- 
firmed by  analysis  of  the  differential  equation  for  various  functional 
forms  of  the  isotherm,  including  cases  where  it  was  not  single-valued 
(45,  46,70). 

As  would  be  expected,  the  assumption  of  instantaneous  equili- 
brium is  not  valid  for  all  situations.  Exact  solutions  of  the 
differential  equation  incorporating  two  important  models  of  rever- 
sible chemical  kinetics  have  been  presented  by  Amundson(3).   These 
models  are 

II  =  k^C  -  k2S  (19) 

for  linear  adsorption  kinetics,  and 

-^^-  =  k,C(S    -  S)  -  k„S  (20) 

9t      1   max         2 

for  "Langmuir"  kinetics,  where  Sj^.^^  ^^  ^^'^^   saturation  capacity. 

Numerical  procedures  have  increased  the  range  of  possible  theoretical 

models  that  can  be  used. 

Except  for  a  few  cases  (24,  69),  less  effort  has  been  devoted 

to  studying  diffusion  kinetics.  In  soils,  however,  where  particles 
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exist  as  aggregates  and  dead  end  pores  are  present,  diffusion  kine- 
tics may  be  important  (69). 

Investigation  of  the  simultaneous  transport  of  a  number  of 
solutes  which  influence  the  mass  transfer  processes  of  each  other 
have  not  been  attempted.   In  addition,  comprehensive  reviews  have 
not  been  made  of  the  various  equilibrium  and  kinetic  models  commonly 
used  by  investigators.   The  equations  derived  from  theory  reflect 
the  effects  of  only  those  interactions  for  which  the  models  have 
been  developed.   Different  theories  may  well  attribute  experimental 
results  to  quite  different  causes.   This  constitutes  the  deficiency 
in  maay  studies  since  it  is  often  extremely  difficult  if  not  impossible 
to  independently  measure  the  theoretical  parameters  introduced  in 
sorption  kinetic  models. 

Equilibria  and  Kinetics  of  Inorganic 
Ion  Exchange  Adsorption 

Inorganic  ion  exchange  was  discovered  more  than  a  century  ago. 
For  soils,  ion-exchange  properties  were  traced  to  their  alumino- 
silicate  fractions  (73,74).   For  many  silicates,  only  ions  from  the 
exposed  surface  layers  are  in  a  position  to  exchange.   Typically 
grinding  increases  the  exchange  capacity  of  these  minerals  (41). 
For  others,  especially  the  zeolites,  the  lattice  structure  is  open, 
permeated  by  waterfilled  channels,  and  possess  internally  accessible 
exchange  sites.   Industrial  applications  of  these  and  other  synthetic 
aluminosilicatcs  were  limited  by  their  low  exchange  capacity  and 
instability  under  acid  or  alkaline  conditions.   The  recent  intro- 
duction of  synthetic  organic  ion  exchangers  with  superior  properties 


led  to  widespread  indastrial  .aid  laboratory  exploitation  of  ion 
exchange  and  a  resultant  surge  in  research  on  ion  exchange  equilibria 
and  kinetics.   With  these  exchangers,  it  was  possible  for  the  first 
time  to  vary  their  properties  systematically.   Much  of  the  theore- 
tical advances  made  were  due  to  this  fact. 

Ion  exchange  equilibria  and  kinetics  are  fundamental  to  an 
understanding  of  ion  transport  and  have  been  studied  intensively. 
As  yet,  no  single  comprehensive  theory  exists  to  explain  all  of  the 
results  involving  ion  exchange.   Certain  characteristics  are  normally 
common  to  all  ion-exchange  reactions.   The  reaction  is  usually  revers- 
ible and  always  involves  an  equivalent  exchange  of  ions.   In  addition, 
the  exchanger  preferentially  adsorbs  one  ion  over  another,  a  property 
appropriately  termed  "selectivity." 

Electrical  double-layer  theory  (71)  advanced  by  Helinholtz  and 
modified  by  Couy  and  Stern  as  an  explanation  of  the  electrokinctic 
properties  of  colloids  has  been  utilized  to  explain  the  phenomena 
associated  with  ion  exchange  by  soil  clays  (7,22,37).   Tliis  tlieo/y 
adequately  accounts  for  the  preferential  sorption  of  ions  of  higher 
valence  over  those  of  lower  valence,  but  fails  to  explain  observed 
exchanger  selectivity  among  ions  of  equal  valence.   Bolt  (8)  attempted 
with  limited  success,  to  extend  the  double  layer  theory  to  account  for 
such  behavior. 

The  overall  ion  exchange  process  can  formally  be  represented  as 
a  reversible  chemical  reaction.   The  exchange  of  ions  1  and  J   of 
valence  +  z..  and  +  z.  respectively  can  be  v^ritten  as 
z.  Il^i|  +  z.  II^Jl   -  z.  l|'i|  +  z.  Jl^Jl 


17 


where  the  bar  signifies  the  ion  associated  with  tlie  exchanger.   As 
dainonstrated  by  Kerr  (A2),  the  ion  in  combination  with  the  exchanger 
could  not  be  considered  as  a  precipitate  with  unit  activity.   He 
assumed  that  the  combination  behaved  as  a  solid  solution,  a  view  that 
was  subsequently  supported  (5,11,72),  and  has  gained  wide  acceptance 
by  investigators  studying  ion  exchange. 

A  useful  quantity,  the  selectivity  coefficient,  characterizing 
the  relative  preference  of  the  exchanger  for  the  ions  I  and  J  was 
obtained  by  applying  the  law  of  mass  action  without  activity  correc- 
tions.  Thus 

lA     =  _J^ J (21) 


'         ^J^i|  .  c.|-j| 


J 

where  C  represents  molar  concentrations" 

Except  when  z-  =  z.,  the  numerical  value  of  K  depends  on  the 

choice  of  concentration  units.   The  total  equivalent  concentration 

Co  =  z   C  +  z   C   of  ions  I  and  J  in  the  solution  and  the  corre- 
i   i    J   J 

spending  quantity  Qo  =  z .  C.  +  z   CT   associated  with  the  exchanger 
must  remain  unchanged  throughout  the  reaction.   Defining  X  ^  z.C./C^ 
and  X.  =  z.C./Q   as  the  "equivalent  ionic  fractions"  of  the  ion  I 
in  solution  and  exchanger  respectively,  and  similar  quantities  for 
the  ion  J,  yields 

X.\''i\    .   X.|-i|      °  (22) 

The  selectivity  coefficient  as  defined  is  thus  dependent  on  the  total 
equivalent  concentration  of  the  exchanging  ions  in  solution  and  the 
capacity  of  the  exchanger.   In  addition,  as  sho\Nm  by  Bonner  and  Bonner 
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ar.d  Payne  (9,10),  it  is  also  a  function  of  the  extent  of  exchange  and 
assumes  its  highest  value  when  the  exchanger  is  completely  saturated 
with  ion  J.   The  selectivity  coefficient  is  thus  not  a  particularly 
useful  quantity  for  predicting  iou  exchange  equilibria. 

A  thermodynanically  rigorous  application  of  the  mass  action  law 
results  in  an  expression  for  the  thermodynamic  equilibrium  constant. 
This  expression  is 

a,|=j|  .  aj|H|  (23) 

where  "a"  represents  ion  activity.   Determination  of  K"  requires 
quantitative  values  for  the  activity  coefficients  of  the  exchanging 
ions  in  solution  and  in  combination  vjith  the  exchanger.   This  con- 
stitutes the  difficulty  in  using  the  classical  thermodynamic  api)roacli. 

II 
The  Debye-Huckel  theory  and  its  extensions  to  predict  activities  of 

electrolytes  in  solutions  are  applicable  for  concentration  ranges 
below  those  normally  associated  with  the  solid  solution  concept. 
This  difficulty  is  not  insurmountable,  hox/ever,  because  methods  based 
on  classical  thermodynamics  are  available  to  determine  the  activity 
coefficients  of  ions  in  combination  v/ith  the  exchanger   and  the 
equilibrium  constant  (17,26).   This  treatment  is  of  more  theoretical 
than  practical  interest  because  it  requires  numerous  measurements 
before  the  equilibrium  constant  can  be  determined  (21,30). 

Ion  exchange  as  a  Donnan  system  was  introduced  by  Mattson  and 
Larsson  (A7)  ,  and  probably  represents  the  most  po\;erful  concept  in 
explaining  ion  exchange.   It  v;as  mainly  through  the  work  of  Bauman 
(3),  Boyd  (11)  and  Glueckauf  (27)  that  this  important  concept  has 
gained  general  acceptance.   The  quantitative  principle  involved  in 
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tlic  Donnan  approach  is  essentially  a  generalization  of  the  double 
layer  concept,  which  is  of  universal  occurrence  v/henever  electrical 
charge  is  confined  within  a  definite  region  of  space.   The  Donnan 
approach  imposes  the  thermodynamic  condition  tliat  all  other  movable 
charges  must  adjust  themselves  accordingly  to  produce  a  niininum  in 
the  free  energy  of  the  system  at  equilibrium.   The  Uonnan  treatment 
thus  represents  a  partial  fusion  of  the  electrostatic  and  thermo- 
dynamic aspects  of  the  ion-exchange  phenomenon. 

Another  approach  likens  the  ion-exchange  process  to  the  phys- 
ical adsorption  of  gases.   Both  the  Langmuir  and  Freundlich  equa- 
tions have  been  used  to  describe  data  on  ion  exchange  equilibria 
(11,3/).   There  is  a  general  concensus  that  ion  excliange  processes 
involve  strong,  long  range  electrostatic  forces  which  were  not  con- 
sidered in  the  conceptual  development  of  the  Langmuir  and  Freundlich 
equations.   Only  vague  physical  meanings  can  be  attached  to  the  para- 
meters obtained  from  the  application  of  these  equations,  and  there- 
fore, they  have  not  provided  much  insight  into  the  adsorption  pro- 
cesses involved. 

A  novel  and  interesting  simplified  statistical  approach  has 
been  used  by  Jenny  (38)  and  Davis  (18)  to  derive  general  theoretical 
equations  for  ion-exchange  equilibrium.   Although  the  results  appear 
interesting,  these  concepts  have  not  been  used  widely. 

Although  ion  exchange  can  be  represented  formally  as  a  chemical 
reaction,  the  physical  processes  involved  have  little  in  comraon  with 
true  chemical  reactions.   Evidence  for  this  has  appeared  in  equilibrium 
studies  (42,72)  where  the  ion  in  combination  with  the  exchanger  was 
treated  as  a  solid  solution.   In  addition,  standard  enthalpy  changes 
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for  ion-exchange  reactions  are  often  less  than  two  kilocalories  per 
ruolc  whicli  is  typical  of  the  orientation  energies  involved  in  dipole- 
dipole  interactions.   Such  evidence  indicates  that  ion  exchange  is 
essentially  a  statistical  redistribution  of  the  exchanging  ions  be- 
tween the  exchanger  and  the  solution. 

Further  support  for  this  concept  was  provided  by  the  pioneering 
iou-oxchange  kinetic  studies  of  Boyd,  Adarason  and  Meyers  (12).   They 
obtained  close  agreement  between  their  experimental  data  on  exchange 
kinetics  and  predictions  from  theoretical  equations  based  on  the  con- 
cept of  ion  exchange  as  a  diffusion  process.   They  hypothesized  that 
either  diffusion  through  a  stagnant  film  around  the  particle  or 
diffusion  into  the  particle  were  rate  controlling.   Integrated  rate 
equations  for  both  cases  were  obtained  by  application  of  Ficks'  laws 
for  constant  diffusion  coefficients.   The  hypothetical  stagnant  film 
was  assumed  to  have  a  finite  thickness  and  was  regarded  as  similar  to 
the  'Nernst'  film  encountered  in  reactions  at  electrode  surfaces,   Tliey 
found  that  either  one  or  both  mechanisms  were  rate-controlling  and 
depended  upon  the  experimental  conditions.   Their  findings  were 
immediately  confirmed  by  other  investigators  (32,4A,58).   These 
studies  showed  that,  in  general,  particle  diffusion  kinetics  were 
favored  for  solution  concentrations  greater  than  10~^  H, efficient 
mixing  (which  reduces  film  thickness)  ,  large  particle  size  and  lov>7  diffu- 
sion coefficients  for  ions  in  the  exchanger.   Opposite  conditions  were 
conducive  to  film-diffusion  kinetics.   Ileasurements  of  self-diffusion 
coef f icietits  of  cations  in  syntlietic  organic  exchangers  were  made  by 
Koyd  and  Soldano  (13).   In  gcnieral,  these  were  an  order  of  magnitude 
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less  than  the  corrGspoiiding  values  in  solution  and  decreased  with 
increasing  valence. 

A  significant  improvement  in  the  analysis  of  the  two  step 
diffusion  concept  was  introduced  by  Helfferich  and  his  coworkers 
(34,36,56).   They  observed  that  Pick's  first  lav/  did  not  take  into 
account  electrokinetic  forces  involved  in  the  interdif fusion  of  two 
charged  species.   They  introduced  the  raore  appropriate  Nernst-Planck 
flux  equation,  which  contained  an  additional  electrical  transference 
term.   This  equation  is 

^     RT  ^^   8x  (24) 

where  J*  is  flux  of  any  charged  species,  F  is  the  Faraday  constant 
and  ^    is  electrical  potential.   Equation  (24)  was  derived  for  diffu- 
sion of  charged  particles  in  an  electrical  field  assuning  ideal 
systems  and  is  widely  used  in  the  analysis  of  electrochemical  reactions. 
The  requirement  that  electroneutrality  be  maintained  at  all  points  in 
the  system  implies  a  rigid  coupling  of  the  concentrations  and  fluxes 
of  the  exchanging  ions.   Thus,  for  two  ions  I  and  J,  electroneutrality 
requires 

|z^|  C.   +  |z.|  C.   =   constant  (25a) 

and 

Z.J*   +  z  f":        =  0.  (25b) 

11       3    2 

These  conditions  allow  the  derivation  of  a  single  flux  equation 
for  either  I  or  J.   This  equation  is  given  (3  6)  by 
J . .  ^    _  DiPjC^i^c.  -f  Zj2c.)   3C, 

Zi2c.D,  -f  Z.2C.D.     ^'^  (26) 

An  immediate  observation  from  this  equation  is  that  the  interdif fu- 
sion flux  is  dependent  on  the  relative  concentrations  of  the  inter- 
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diffusing  ions.   With  vanishing  concentration  of  either  I  or  J 
the  interdif fusion  flux  is  controlled  by  the  diffusion  coefficient 
of  the  ion  in  the  minority.   Further  analysis  of  the  kinetic  be- 
havior of  ion-exchange  reactions  by  these  investigators  for  both 
film  and  particle  diffusion  (36)  led  to  the  conclusion  that  the 
rates  of  forward  and  reverse  exchange  reactions  were   not  equal. 
This  conclusion  vjas  confirmed  experimentally  (34,35). 

Incorporation  of  the  simple  diffusion-kinetics  model  based  on 
Pick's  ]aw  into  the  differential  mass-transport  equations  pro- 
duces a  mathematical  problem  of  extreme  complexity  (40).   As  a 
result  simpler  equp.tions  based  on  the  linear  diffusion  concept 
introduced  by  Glueckauf  (29)  have  been  utilized  (63).   As  pointed 
out  by  Helfferich  in  his  comprehensive  monograph  on  ion  exchange 
(35),  for  practical  situations,  the  gain  in  accuracy  does  not  warrant 
the  time  and  effort  expended  in  solving  the  problem  of  greater  com- 
plexity. 


CRr\PTi:R  2 

SYNTHESIS  OF  PERTINENT  THEORY  ON 
INORGANIC  ION  EXCHANGE  AND  TR/iNSPORT  PROCESSES 

Thernodynamic  Conceptualization  of 
Ion-Exchange  Processes 

x\  system  is  defined  thernodynamically  as  a  body  or  group  of  in- 
teracting bodies  intended  for  separate  study.   Any  physically  homo- 
geneous body  or  set  of  identical  homogeneous  bodies  is  called  a  phase. 
Phases  are  either  pure  or  mixed,  depending  on  whether  they  consist  of  a 
single  or  several  chemically  individual  species. 

Systems  are  either  homogeneous  or  heterogeneous  depending  on 
vliether  they  consist  of  a  single  or  several  phases.   The  existence  of 
physical  boundaries  (interfaces)   and  interphase  regions  are  necessary 
features  of  polyphase  systems.   A  solution/exchanger  system  can 
therefore  be  conceived  thermodynamically  as  a  heterogeneous  system  con- 
sisting of  two  mixed  phases. 

A  heterogeneous  system  may  exist  either  in  an  equilibrium  or  non- 
equilibrium  state.   In  the  former  state  all  thermodynamic  state 
variables  remain  constant  with  time.   If  the  system  is  non-equilibrium, 
spontaneous  phase  interactions  occur  resulting  in  the  establishment  of 
an  equilibrium  state,  characterized  by  definite  compositions  of  all  the 
phases.   Interactions  that  do  not  involve  the  production  of  new  phases 
or  new  chemical  compounds  result  in  material  or  energy  exchanges 
across  the  interfaces.   Such  interactions  are  considered  as  physical 
sorption  processes   and  involve  atomic  and  molecular  interaction 
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energies  distinct  from  tlior.e  involved  in  cheraica]  bonds.   The  term 
adsorption  refers  to  physical  sorption  in  which  the  species  trans- 
ferred becomes  either  concentrated  at  the  interface  or  distributed  in 
the  bulk  of  the  phase.   If  instead  of  being  transferred,  tlie  species 
is  displaced  by  interactive  forces  back  into  the  same  phase  it  is 
termed  negative  adsorption  or  loosely  as  exclusion. 

A  unique  property  of  exchanger  phases  in  the  presence  of  fixed 
electrical  charge  sites,  which  may  be  either  restricted  to  the  ex- 
changer surface   or  distributed  throughout  its  bulk  volume.   The 
quantity  of  fixed  charges,  conveniently  expressed  as  equivalents, 
defines  the  absolute  capacity  of  the  exchanger.   This  property, 
coupled  with  the  restriction  that  electroneutrality  be  satisfied  at 
all  points  in  either  phase  of  the  system,  forms  the  basis  for  the 
metathetical  sorption  phenomena  in  solution/exchanger  systems „   Uere  it 
not  for  its  fixed  charges,  the  exchanger  would  lose  its  identity  as  a 
distinct  phase  in  the  system  at  equilibrium.   Ion  exchange  processes 
occurring  during  equilibration  of  a  solution/exchanger  system  can  be 
considered  as  physical  adsorption,  if  it  is  hypothesised  that  no  in- 
teractions occur  involving  the  formation  of  chemical  bonds  and  produc- 
tion of  new  chemical  species  or  phases  in  the  system. 

Disposition  of  Charged  Species  in 
Solution/Exchanger  Systems  at  Equilibrium 

For  charged  species  in  a  heterogeneous  system,  a  necessary  con- 
dition at  equilibrium  is  equality  of  the  electrochemical  potential  of 
each  species  in  the  various  phases.   The  electrochemical  potential  n 
for  a  species  in  a  phase  is  defined  (31)  by 

n  =  li  +  zF<))  (27) 
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in  wliicli  p  is  the  chemical  potential,  z  the  valence,  F  the  Faraday 
constant  and  i|)  t)ie  inner  potential  of  the  phase.   The  electrochemical 
potential  can  be  conceived  as  the  sum  of  the  reversible  chemical  and 
electrical  work  required  to  transfer  a  charged  particle  from  infinity 
to  any  point  in  tlie  interior  of  the  phaseo   If  the  particle  is  un- 
charged no  electrical  work  is  involved  and   n  =-\io 

Consider  two  ions  I,  J  of  valence  z.,  z.  in  the  solution/ 
exchanger  system  at  equilibrium.   Then,  r\ .    =  n .  and  n.  ~  n.>  ^^here  the 
bar  signifies  the  exchanger  phase.   The  chemical  potential  for  a 
species  in  a  phase  is  given  by 

y  =  y°  +  RT  iln  a  (28) 

where  y"  is  the  chemical  potential  in  an  arbitrary  reference  state,  R 
is  the  gas  constant,  T  is  absolute  temperature  and  "a"  is  activity. 

Substituting  equation  (28)  into  equation  (27)  and  equating  the 
electrochemical  potentials  of  I  and  J  gives 

y?   +  RT   Jin  a.    +  z.    F  "^  =   \i°  +  RT   Zn   a .   +   z  .    F    <>  (29a) 

y?   +  RT   £n  a.   +  z.    F  ^  =  y!  +  RT   iln  a.   +   z.    F   0  (29b) 

The  difference  in  the  inner  potential  of  the  phases  ^  -  (J<,  at 
equilibrium,  is  invariable,  and  \i°   =   V° ,    y?  =  P?.   Thus 

-    ,    RT  „   ai    RT  „   ai  ,„„  . 

*  -  *  =  TT  '"  i^  =  7J  ^"  f"         ^29c) 

and 


Multiplying  by  z.  z.  equation  (29d)  becomes 


(29d) 


z.  ;,n  ^  =  z.  £n  ^  (29e) 

J    a.     .    a. 


•whenci?. 
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(£i)2j  =  (^fi  (29f) 

a.       a . 
1       3 


Similar  reanoning  for  an  ion  X  of  valence  ~z^   yields 

(30) 

XI  X  J 

Insight  into  the  usefulness  of  the  above  relations  can  be  ob- 
tained by  considering  some  specific  cases. 

Consider  an  exchanger  with  fixed  negative  charges   and  absolute 

capacity  Qo  satisfied  by  J  ions  in  equilibrium  with  a  solution  of  an 

electrolyte  J  X   o  Electroneutrality  requires 
^x  Zj 

z.C.  =  Qn  +  z  C   and   z.C.  =  z  C  ,         (31a) 

JJ°      XX  JJ      XX* 

where  C  represents  the  molar  concentration.   From  above   equation 
(29c)  shov7S  that 

*  -  ♦  ■=  S  ""  f^ 

J      3 
The  potential  difference  ^  -  ({>  is  the  equilibrium  Donnan  potential 
across  the  interphase,  and  increases  ai>/ay  from  the  solution/exchanger 
interface.   The  interphase  functions  as  a  Donnan  membrane  in  a  ther- 
modynamic sense  because  it  is  impermeable  to  the  fixed  exchanger 
charges.   This  macroscopic  potential  is  immediately  established  and  is 
tb.e  force  preventing  the  net  transfer  of  J  ions  out  of  the  exchanger 
and  of  X  ions  into  the  exchanger  despite  concentration  differences  that 
exist  between  the  two  phases.   Increasing  concentrations  of  J  ions  in 
solution  causes  a  lowering  of  the  Donnan  potential  while  ir.croasing 
the  exclianger  capacity  \i7ould  result  in  larger  potentials.   It  is  ob- 
vious that   lim  (^--(J))  -■  0.   If  a  dilute  solution  J-,  Xy .    is  used  for 

aj-.clj  _    _  _         '""      '^ 

equilibration  tlion  a.  -  C.   and   a.  =  Y.C.,  where  y   denotes  the  activity 
'  ^  2  3     3  J 

coef ficiento   Irom  equation  (31a)  it  follov;s  that 
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Qo  +  zy  Cx  ^  So  /3^,^N 

z.       z. 
2  J 


From  above  equation  (30)  is 


(f^)^3    =  (^)^x 
a        a . 
X        J 

Substituting  C.  for  a.  and  ^3^     for  a   gives 
J      J      Zj        3 

(ix)^j   .   (|^)^x  •    (31c) 

X         '  ^o 

Since  y.Q  is  much  larger  than  C.z.,  the  concentration  of  X  in  the 
exchanger  is  lower  than  in  solutiouc   This  effect  is  greater  the  more 
dilute  the  solution  of  J^  X„ . ,  and  is  enhanced  with  an  increasing 

Zjj  Zj» 

valence  of  the  ion  X„   It  is  commonly  called  "Donnan  exclusion." 

If  a  salt  I7  X„.  is  now  introduced  into  the  system  described 
•^x  ^x 

above,  exchange  of  I  and  J  occurs^   At  equilibrium  equation  (29f)  shows 
that 

a.       a^ 


Rearranging  gives 


(ii)^j  ^  (ai)^J 


(a^)'i    (a^)' 


1  (32) 


If  I  is  introduced  in  a  trace  quantity,  then  a.  is  approximately 

constant.   With  this  condition,  application  of  Le  Chatelier's  Principle 

to  equation  (32)  shows  that  if  the  concentration  of  J  is  increased  in 

the  equilibrium  solution  then  a.  decreases  and  vice-versa.   Thus,  for 

dilute  solutions  of  I2  X^ .  the  presence  of  J  ions  in  the  solution 

phase  will  lov/er  the  selectivity  of  the  exchanger  for  I  ions>.   This 

effect  will  be  greater  the  larger  the  difference  between  z.  and  z  o 

From  the  general  relation  given  by  equation  (29f) 

(f  )^j  =  (|i)^i 
a.       a^ 
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Dividing   thi-ougliout   by    (=>^)    ■"-   results    in 


(33a) 
-i  -j 

Introducing  activity  coefficients  gives 

The  terra  on  the  left  side  of  equation  (33b)  is  the  definition 
of  the  "selectivity  coefficient"  obtained  by  application  of  the  nass 
action  law  v/ithout  activity  corrections  for  the  exchange  of  I  and  J„ 
Since  the  activity  coefficients  are  functions  of  species  concentration, 
the  selectivity  coefficient  depends  on  the  experimental  conditions. 

Physical  Basis  for  Exchanger  "Selectivity" 
The  power  of  thermodynamics  is  its  ability  to  produce  general 
relations  among  system  variables  without  detailed  knowledge  of  the 
specific  physical  forces  involved  in  the  phase  interactions.   These 
are  coiicealed  in  the  thermodynamic  activities  of  the  system  components. 
The  equilibrium  composition  of  the  phases  is  governed  by  the  require- 
ment that  the  free  energy  of  the  system  be  minimized „   A  charged 
species  in  a  non-equilibrium  solution/exchanger  system  can  lower  its 
free  energy  by  interactions  with  the  exchanger,  solvent,  and  other  ions 
in  Che  system.   Quantitative  theories  do  not  exist  to  predict  exactly 
the  effect   of  such  interactions,  generally  termed  solvation  processes, 
on  the  thermodynamic  properties  of  the  phases.   The  sum  total  of  all 
interactions  reveals  itself  in  the  observed  selectivity  property  of 
tlie  exchangerc   In  effect,  selectivity  is  a  measure  of  the  relative  in- 
teractive effect  of  two  exchanging  ions  on  the  thermodynamic  properties 
of  the  exchanger  phase. 
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In  an  aqueous  solution/exchanger  system,  electroneutrality  re- 
quircrf  that  the  fixed  charge  on  the  exchanger  be  satisfied,  at  all 
tiraes,  by  an  equivalent  quantity  of  charges  of  opposite  sign.   This 
fact  deterraines  the  minimum  equivalent  concentration  of  the  intersti- 
tial solution  of  porous  exchangers.   If  Q  represents  the  absolute 
capacity  in  equivalents  per  gram,  p  the  dry  bulk  density,  and  f  the 
fractional  internal  porosity  of  the  exchanger  particles,  then  pQ  /f 
equivalents/cm  is  the  concentration  of  the  interstitial  solutiono   An 
exchanger  with  Q  =  100  meq/lOOg,  p  =  1  g/cm  ,  and  f  =  0.5  gives  a 
concentration  of  2  N  for  the  interstitial  solution.   Similar  concen- 
trations would  occur  in  the  interphase  regions  of  exchangers  with  sur- 
face charge  sites ^   In  high  capacity  synthetic  organic  exchangers, 
values  as  high  as  10  N^  are  encountered o   Such  concentration  ranges  are 
beyond  the  scope  of  the  Debye-Hiickel   theory  for  obtaining  the  thermo- 
dynamic effect  of  ion-ion  interactions c 

Ion-solvent  effects  can  be  evaluated  by  the  Born  theory.   This 
theory  gives  the  solvation  free-energy  per  mole  of  ions  in  solution 


^       2(r  +  0.85)  ^^        c^  ^      ^ 


e  the  electron  charge,  e  the  dielectric  constant,  and  0o85  A  an  em- 
pirical correction  factorc   The  change  in  free  energy  in  transferring 
one  mole  of  I  ions  from  solution  to  an  exchanger  in  the  J  form  due  to 

ion-solvent  interactions  and  electrical  work  would  be 
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Since  tha  exchanger  solution  is  nore  concentrated  c  <  e  and  the 
solvent  interaction  term  is  positive „   Thus  the  transfer  does  not 
occur  spontaneously  unless  the  second  term,  which  is  negative,  has  a 
greater  absolute  valuCo 

The  transfer  of  I  into  the  exchanger  must  be  accompanied  by  an 
equivalent  transfer  of  z./z.  moles  of  J  out  of  the  exchanger,,   For 
this  process 

'''  '  ul'^'liy  <f  -  i>  *  ^/'*  -  *>  "^'=' 

In  this  case  the  leading  term  is  negative  and  the  electrical  work  term 
is  positive^   The  total  free  energy  change  for  the  metathetical  reac- 
tion is 

AG  =  AG  +  AGo  (34d) 

Wien  z.    =  z.,  AG  is  negative  if  r.  <  r .  c,   Thus,  considering  only  ion- 
solvent  interactions,  an  exchanger  with  fixed  charges  neutralized 
by  J  has  preference  among  ions  of  equal  valence  which  would  increase 
in  order  of  increasing  crystal  radius <,   In  part,  ion-solvent  interac- 
tions serve  to  explain  selectivity  among  ions  of  equal  valencec   V.Tien 
z.  4   z.,  the  electrical  work  term  dominates  the  solvation  terra  and 
thus  regardless  of  radius,  the  exchanger  prefers  the  ion  of  higher 
valencGc 

The  concept  of  exchanger  selectivity  as  competitive  solvation 
cannot  be  developed  further  because  a  complete  understanding  of  these 
processes  is  far  from  being  realized,   Hov;ever,  the  above  does 
illustrate  the  complexity  of  the  physical  interactions  involved  in 
solution/exchanger  systems , 
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Concepts  on  Inorganic  Ion  Transport 
Durinf;  t'iscible  Displacenent 

Established  concepts  on  ion-exchange  equilibria  and  kinetics 
allow  some  qualitative  deductions  regarding  the  transport  behavior 
of  inorganic  ions  during  miscible  displacements. 

If  Instantaneous  exchange  equilibriura  is  assuned,  the  transport 
beliavior  of  an  ion  depends  on  the  shape  and  characteristics  of  the 
adsorption  isothem.   These  would  rest  heavily  upon  the  selectivitj/ 
properties  of  the  exchanger.   Over  a  snail  range  of  very  dilute  con- 
centrations of  the  equilibrating  solution  phase,  it  can  be  expected 
that  the  isotherms  \;ould  be  linear.   In  such  ranges,  for  an  exchanger 
saturated  with  J  ions,  the  presence  of  J  ions  in  the  equilibrating 
solution  would  influence  the  sorption  of  another  counterion  in  the 
system.   As  a  result,  variations  in  the  concentration  of  J  ions  in  a 
displacing  solution  of  I  ions  would  produce  variations  in  the  elution 
tines  for  the  ion  I. 

The  two-step  diffusion  concept  of  exchange  kinetics  suggests  tliat 
for  snail  diameter  exchanger  particles  and  trace  concentrations  of  I 
in  the  displacing  solution,  film-diffusion  kinetics  nay  control  mass 
transfer  of  I.   Since  the  film  is  considered  as  a  stagnant  hydro- 
dynamic  boundary  layer  around  the  exchanger  particle,  its  thickness 
would  be  inversely  influenced  by  flo\/  velocity  of  the  solution.   Also, 
as  discussed  in  the  foregoing  review  [equation  (26)],  mass  transfer  of 
I  by  diffusion  would  be  influenced  by  the  presence  of  J  ions.   The 
existence  of  such  rate-controlled  diffusion  processes  would  reflect 
in  deviations  from  the  predicted  equilibrium  shapes  of  elution  curves 
for  I. 


CHAPTER  3 

SOLUTION  OF  THE  COOTECTIVE-DISPERSION 
EQUATION  FOR  A  "PULSE  INPUT"  AND  ITS  EXTENSION 
FOR  A  SOLUTE  FOLLOl/ING  A  LINEAR  ISOTHEMI 

Although  solutions  to  the  convective-dispersion  equation  arc 
often  quoted  and  used,  detailed  derivations  of  these  solutions  are  not 
usually  presented.   Details  of  an  asj^ptotic  solution  to  the  convective- 
dispersion  equation  for  "pulse  input"  boundary  conditions   are  given 
below  for  a  non-reactive  and  a  reactive  solute  (linear  adsorption 
isotherm) . 

The  raathematical  formulation  of  the  problem  is  given  by  the 

equation 

2 

-1^-  =  D^  -  V  1^  for  0  <  X  <  CO  and  t  >  0       (35) 
9t    dx'^  9x 

with  initial  and  boundary  conditions 

C(x,o)  =0,  for  0  <  X  <  <»  (36a) 

C(0,t)  =  Co[U(t)  -  U(t  -  t^)],  for  t  >  0      (36b) 

C(x,t)  =  0,  for  t  >  0  (36c) 

X  ->  CO 

v/aere  U(t)  is  the  Hcaviside  unit  function  and  t,  is  the  time-width  of 
the  pulse., 

Let  tlie  Laplace  transform  of  C(x,t)  be  denoted  by  u(x,s)c   The 
above  equation  and  auxiliary  conditions  under  the  Laplace  transforma- 
tion become 

su(x,s)  -  C(x,0)  =  D^-  v^.  (37) 
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u(0,s)  =  -^(1  -  e-^^1).  (38a) 


u(x,s)  =  0.  (38b) 


The  characteristic  equation  of  the  ordinary  differential  equation  is 


Dm  -  vm  -  s  =  Oo  (39a) 


with  roots  v  +  /\'-  +   4Ds. 
2D 


Denoting   /v-  +  ADs  hy  R,  the  general  solution  is 
2D 

u(x,s)  =  C^e<IE  -^  ^)^  +  C^e^i^  "  ^^\  (39b) 
Applying  the  above  initial  and  boundary  conditions 

u(0,s)  =  C-^  +  C2  =  -^(1  -  e~^h)„  (39c) 

u(x,s)  =  0  =:^  C^  =  Oo  (39d) 

X  ->  00 


u(x,s)  =  -^(1  -  e~'"-l)e'^  ~  ^^^^„  (39e) 

Co^(^-R)x  Coe"^^l   e^^-^)^.  (39f) 

s  s 

C  e'^  "    e:^  -  C  e^  «  e"^^l  e"^^  (39g) 


Consulting  a  table  of  Laplace  transforms  (60),  transform  pair  //3i,2-80  is 

listed  as 

1    /  ^  T,2.1/2     -ab  ,  ,_ 

1  -a(s  +  b  )  ^    ^  e       c  /  3      i,^l/2x 
-^         '< *  -^—     erfc(- ^  -  bt    ) 

vaJid   for   real   s   >0  » 


Since 


-   „-,/^+  4Ds,x  -  Tm//,^/^2 


V  4D 


(40) 


,    the   first      terra  inverts   to 

D  2/W 


vx     vx              1/2     vx              1/2 
, )  +  e^U  erfc( +  

2v^        2/T  2v^        2/^" 


^  e2D   {e      2D  erfc(-^  -  ^^  )   +  a'^  erfc(-^  +  ^^£-^)}(41a) 


=  ^  {erfc(i^^^)   +  0-15-  erfc(^^-i^)}  (Alb) 

2/Dt  2/Dt 

The  second  term  is  inverted  by  noting  that  it  is  equal  to  the  first 

—  St 

term  x  e    >.   By  the  shifting  property  of  the  Laplace  transforms,  if 
L{f(t)}  =  f(s),  and  g(t)  =  f(t  -  ti).   U(t  -  t^^)  then,  L{g(t)} 

— St  1 

=  e     f(s)c   The  second  term  inverts  to 

f  {erfc["-"(^-  ^^]  +  e^  erfc[^-^  "  ^'  '   '^l]}    ^    U(t  -  t^)  (42) 
2/D(t  -  ti)  2v^D(t  -  ti) 

The  second  term  in  both  inverse  transforms  is  small  except  near  the 

inlet  wliere  x  =  0  and  for  small  values  of  t  unless  D  is  large.   As  a 

result  it  can  be  ignored  without  introducing  a  serious  erroro   The 

solution  reduces  to 

-^  =   l/2{erfc(^^-^^)   -   erfc    [^  "  "^^   "   ^^^J    U(t  -   tj}  (43a) 


"o  2/  Dt  2/D(t   -    t    ) 

Using   the   identity   erfc    (x)    =   1   -   erf    (x) ,    this   can  be  written  as 

JL  =   i/2{erf[^  -  ^^^   -   '^]    .    U(t  -   t,)   -   orH^^^^^))  (43b) 


2v^(t  -  t^)  -^  i/YF 


It  is  convenient  to  use  a  transformed  variable  G  =  -r-,  which  is 
])hy£ically  equivalent  to  the  number  of  pore  volumes,  for  experiments 
with  a  fixed  value  of  x  =  L.   Setting  x  =  L  and  dividing  the  argument 
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of  the  error  functions  in  equation  (A3b)  top  and  bottom  by  L,  the 
solution  given  by  equation  (A3b)  transfoins  into 


For  a  solute  following  a  linear  isotherm,  S  =  KC,  and  —  =  R: —  o   The 

ot  dt 

differential  equation  becomes 

IP-  =  T>^    9C  _  pK  _8C  ^ 

3t   ^3  "  ^   e   9t  '  ^'^   ^^ 


9C 
Transposing  -r—  and  rearranging 
d  t 


dt        R 


ill  WW 


where  R  =  1  +  -^-  • 
^  V 

For  the  same  initial  and  boundary  conditions  it  is  obvious  that  the 
solution  to  equation  (44b)  can  be  obtained  by  setting  D  =  D/R  and 

V  =  v/R  in  the  solution  given  by  equation  (43b) o   This  yields  after 

vt 
introduction  of  G  =  —  and  rearrangement 

S.  .  i/2{erf[IS-=-l±JlI^]   u(e  -  e,) 


2/D(e  -  ei)R 

.   erf[<^-/)    •    ^-]} 
2/     DOR 


(45) 


Equation  (45)  reduces  to  equation  (43c)  when  R  =  1^   The  former 
equation  can  therefore  be  used  to  generate  theoretical  elution  curves 
for  both  reactive  and  non-reactive  solutes,  with  given  values  of  the 
parameters  v,  L,  D,  R  and  6^ .   Equation  (45)  was  obtained  assuming  a 
serai-infinite  space  domain  but  it  can  be  used  to  analyse  displacement 
experiments  in  finite  columns^   It  has  been  shown  (14,  46,  54)  that 
for  large  values  of  v  and  L  no  serious  error  is  involved  in  using  equa- 
tion (45)  for  this  purpose^ 


CHAPTER  4 
EXPERIMENTAL  OBJECTIVES,  IIATERIALS  AND  METHODS 

The  primary  objective  of  this  study  was  to  examine  the  transport 
behavior  of  selected  inorganic  catioi^s  during  steady  saturated  or  un- 
saturated flow  in  a  reactive  porous  medium.   It  was  anticipated  that 
these  results  vv'ould  provide  insight  into  the  nature  of  the  metathetical 
mass  transfer  processes. 

From  the  onset,  it  was  clear  that  soils  were  too  heterogeneous  to 
study  exchange  processes  in  detail.   Therefore,  a  mixture  was  prejjared 
with  a  synthetic  exchange  material  and  sand.   Synthetic  exchangers 
were  utilized  previously  in  miscible  displacement  studies  by  Day  and 
Forsythe  (20).   In  addition  to  providing  the  required  hoiaogeneity,  the 
overall  exchange  capacity  could  be  controlled. 

A  rigid,  analytical  grade,  macroporous,  granular  (50  -  100  mesh), 
organic  exchange  resin  of  a  highly  crosslinked  sulphonated  copolymer 
of  styrene  with  divinylbenzcne  was  purchased  (Bio-Rad  Laboratories). 
This  material  is  both  thermally  and  chemically  stable  and  has  a  cation 
exchange  capacity  of  4.9  meq/g.   In  addition,  special  treatment  during 
polymerisation  results  in  a  low  resistance  to  intra-particle  mass  trans- 
fer by  diffusion.   Preparation  of  synthetic  organic  exchangers  in  gen- 
eral has  been  described  by  Helfferich  (35).   Their  general  physical  prop- 
erties have  been  reviewed  by  Meyers  et  al.  (48)  and  the  special  physical 
and  chemical  properties  of  the  macroporous  resins  have  been  described 
by  Mi]ler  et  al.  (49,50).   Because  of  its  high  excliange  capacity,  it 
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war,  necessary  Lo  dilute  this  resin  with  an  inert  material.   Soil  from 
the  1  m  -  1.3  m  horizon  of  a  Lakeland  sand  (Typic  QuartzipsamruenL) 
was  passed  through  a  nest  of  sieves  and  the  fractions  retained  on  the 
500,  250  and  105-micron  sieves  were  combined  in  the  ratio  of  25  :  50  : 
25,  respectively.   This  material  was  then  treated  with  hydrogen  peroxide 
to  destroy  any  organic  matter  present.   The  resin  was  treated  repeatedly 
with  ]  N^  calcium  acetate  solution  until  no  further  pH  change  was  observed, 
and  then  packed  into  a  plastic  column  and  eluted  with  1  N  Ca(NO  ) 
solution.   This  procedure  was  considered  sufficient  to  Ca-saturate 
the  resin.   The  resin  was  then  dried  and  sieved  and  the  fraction  betvzeen 
200  and  105  microns  combined  with  the  sand  to  yield  a  computed  exchange 
capacity  of  30  -  40  meq/100  g.   This  porous  exchange  material  was  used 
in  all  subsequent  studies. 

Cations  selected  for  investigation  were  Li"^  and  Na  because  of 
their  ease  of  detection  at  low  concentrations  by  flame  spectrophoto- 
metry.  Calcium  v;as  chosen  as  the  common  ion  because  the  Ca-saturated 
exchanger  was  expected  to  exhibit  low  selectivity  for  Li   and  Na 
resulting  in  low  residence  times.   Because  of  this,  longer  columns 
could  be  used.   The  exclusion  of  CI  by  the  exchanger  was  also  studied 
and  a  cursory  investigation  was  made  of  the  transport  behavior  of 
^5ca2+. 

The  columns  to  contain  the  exchange  medium  were  prepared  from  a 
single  length  of  rubber-cast  plexiglass  tubing  with  an  internal  dia- 
meter of  5.0  cm.   The  columns  were  designed  to  provide  unsaturated 
flow  conditions.   To  achieve  unsaturation,  3-mm  holes  were  drilled 
in  the  walls  of  the  column,  and  small  pieces  of  wire  gauze  (less  than 
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100--n:".cron  mesh)  were  placed  on  the  inside  to  retain  the  material  in 
the  column.   Porous  ead-plates  with  an  air-entry  pressure  of  30  -  40 

cm  o C  water  wore  used  to  retain  the  material  tightly  in  the  columns. 

3 

The  dead  volume  of  the  end  plates  did  not  exceed  7  cm  .   The  holes  in 

the  column  were  plugged  during  saturated  flow  studies.   For  unsaturated 

flow  studies  the  holes  were  not  plugged  and  the  column  was  sealed  into 

a  10-cm  diameter  column.   The  outer  column  was  pressurized  to  achieve 

the  desired  unsaturated  soil-water  potential  and  the  pressure  was  kept 

constant  with  a  bubble  tower. 

The  solid  matrix  material  was  packed  into  the  column  under  water 

to  insure  complete  water  saturation.   An  adjustable  peristaltic  pump 

v;as  used  to  deliver  solution  at  predetermined  rates  to  the  columns. 

A  fraction  collector  was  used  to  sample  the  effluent  at  equal  time 

intervals.   A  schematic  of  the  experimental  apparatus  is  given  in 

Figure  1. 

Molar  stock  solutions  of  Ca(NO^)  ,  NaCl ,  LiCl ,  LiNO   and  NaNO 

3  2'      '      '     3         3 

were  prepared.   Tracer  solutions  were  prepared  by  combining  suitable 
portions  of  the  required  stock  solutions  and  making  up  to  volume. 
The  concentration  of  Na  ,  Li   and  Cl~  in  the  tracer  solutions  were 
kept  between  80  -  85  ppm.   Concentrations  of  the  common  ion  (Ca   ) 
used  were  0,05  M,  0.02  M  and  0.005  M.   Tritiated  water  (HTO)  was  added 
Ts  a  non-reactive  tracer  to  evaluate  Viydrodynamic  dispersion.   The 


dissolving  one  gram  of    CaSO,  in  2  liters  oF  0.075  M  and  0.05  M 

Ca(NO  )  . 
3  2 

Exchange  isotherms  for  Na"*"  and  Li"*"  in  0.05  N,  0.02  N  and  0.005  N 
Ca(KO,.)   were  determined  over  a  0  -  200  ppm  concentration  range.   Ten 
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solutions  in  increments  of  20  ppm  were  prepared  by  weighing  out  the 
required  amounts  of  2000  ppm  stock  solutions  and  making  it  up  to  volume 
with  the  appropriate  Ca(NO  )2  solution.   Ten-gram  portions  of  the  dry 
solid  matrix  material,  previously  washed  with  deionized  v/ater,  were 
mixed  with  10-ml  portions  of  each  solution  and  shaken  at  frequent 
intervals  during  a  12-hour  period.   This  time  was  sho^vm  to  be  sufficient 
for  equilibration.   A  sample  of  the  supernatant  was  then  withdrai>m  and 
analyzed.   The  amount  adsorbed  was  calculated  from  the  concentration 
difference  between  the  sample  and  original  solution.   A  similar  tech- 
nique was  used  to  determine  exchange  isotherms  for    Ca    in  0.075  M 
and  0.05  M  Ca (NO  )  . 

The  \)ulse  input'  boundary  condition  was  utilized  in  all  miscible 
displacement  studies.   The  columns  were  leached  with  appropriate 
tracer -free  solutions  of  CaCNO^)^  and  then  a  pulse  of  the  tracer  solu- 
tion containing  the  cation  plus  chloride  and  HTO  was  introduced.   This 
pulse  was  subsequently  eluted  with  the  tracer-free  Ca (NO  )   solution. 
Changeover  of  solutions  was  achieved  in  approximately  2  minutes  during 
which  time  the  outlet  was  sealed,  and  the  front  end-plate  and  delivery 
tubing  were  flushed  and  refilled  with  the  new  solution.   The  total 
aniount  of  solution  introduced  during  a  displacement  was  determined  by 
V7eight  differences  in  the  bottles  containing  the  solutions.   These 
together  with  the  measured  time,  were  used  in  computing  an  average 
pore-v7nter  velocity. 

Analysis  for  Li   and  Na   in  the  effluent  samples  were  made  using 
a  Beckman  model  R  flame  spectrophotometer.   Sets  of  standard  solutions 
for  these  analyses  were  prepared  by  accurate  dilution  of  a  2000  ppm 
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stock  with  the  appropriate  Ca(NO  )   solutions.   Analyses  for  Cl~  were 
made  using  an  "Orion"  specific-ion  electrode.   Activity  of  HTO  and 
^•''Ca^'^  was  determined  by  liquid  scintillation  counting  of  l-ml  samples, 
in  10  ml  of  a  commercial  phosphor  (Aquasol  II). 

The  pore  volume  of  the  columns  was  determined  by  drying  the  material 
in  the  column  at  the  end  of  a  series  of  displacements.   The  dispersion 
coefficients  v/ere  extracted  from  the  breakthrough  data  for  HTO  using 
a  least  squares  curve-fitting  procedure  and  the  simplified  asymptotic 
solution  [equation  (43c)].   These  coefficients  were  then  used  to  gen- 
erate analytical  curves  for  the  reactive  solutes  using  the  sorption 
parameters  from  the  adsorption  isotherms. 


CMPTER  5 

RESUUS  AND  DISCUSSION  OF  STUDIES  ON  EXCHANGE  EiJUILIBRIA 
/lInD  transport  of  Na+,  Li""",  ^^Ca?-^   AND  Cl" 

Exchange  Adsorption  Isotherms 

The  adsorption  isotherms  for  Na"^  and  Li"^  in  0.05  II,  0.02  M  and 

0.005  M  Ca(NO„)„  on  the  exchange  mixture  are  given  in  Figures  2  and  3, 

The  adsorption  isotherms  were  linear  for  both  Na  and  Li  over  the  0  to 

200  ppm  concentration  range.   The  slope  of  the  adsorption  isotherra  in- 

2+ 
creased  as  the  Ca   concentration  in  the  equilibrating  solution  de- 
creased. 

The  data  in  Figures  2  and  3  were  fitted  by  the  least  squares  pro- 
cedure to  the  equation  for  a  linear  isotherm,  S  =  K  C,   The  resulting  K 
values  for  Na"^  in  0.05  M,  0o02  H  and  0.05  M  Ca(N0  )   v;ere  0.1324, 
0.2196  and  0.3770,  respectively.   The  corresponding  values  for  Li  were 

0.0822,  0.1122  and  0.153A.   These  K  values  increase  in  a  non-linear 

2+ 
fashion  with  decreasing  Ca   concentration  in  the  equilibrating  solu- 
tion.  As  sho\>m  in  Figure  A,  an  assumed  general  exponential  relation- 
ship of  the  form  K  -  K  exp  [3  C   ]  fitted  the  data  reasonably  v/ell. 

For  Na  ,  the  values  of  K  and  3  were  0.3883  and  -9.7  and  for  Li   the 
o 

corresponding  values  were  Ool568  and  -5.80.   These  values  provide  use- 

2+ 
ful  quantitative  insight  into  the  damping  effect  of  the  Ca    on  the 

adsorption  of  Na   and  Li  . 

From  the  theoretical  considerations  discussed  previously  [equation 

(32)]  for  ccjuilibriuia  conditions  the-  following  relatinnshi])  holds: 

42 


43 


ouoOLOOLnomoLO 

to  I 

Z2 


Ufy 


A5 


P 

I 


46 


-2 


(46) 


For  the  concentrations  used  in  this  study,  both  a.,  and  a„  are  accessi- 

Na      Ca 

bin  from  the  extended  Debye-Huckel  theory.   This  theory  allows  calcula- 
tion of  the  activity  coefficient  from  the  following  relationship: 

Az^  /I 


-  log  Y  -  Y^r^ 
where  y  is  the  activity  coefficient,  I  = 


(47) 


1  2 

"V  EC.z.  and  is  the  ionic 

2  i  1  X 


strength  of  the  solution,  A  =  0„507  at  20  C,  B  -  0.328  at  20  C  and  a  is 

+       2+ 
the  ion-size  parameter.   The  value  of  a  is  6  for  Li  and  Ca   and  4  for 

Na   (43).   If  only  small  quantities  of  Na  are  adsorbed  in  the  exchanger 

phase  one  can  assume  that  a   is  constant.,   If  it  is  further  assumed 


that  Y,,   remains  constant,  then  C, 
'Na  '       Na 

equation  (46)  to  give 

aNa  ^  Pj:^ 


pSNa 


may  be  substituted  into 


Rearranging 


'Ca  ^^ 


; -—    1  Na 


(48a) 


(48b) 


■'Na 


(48c) 


If  K  is  known,  the  K  value  for  each  adsorption  isotlierm  can  be  computed. 

It  is  recognized  that  since  I  varies  v;ith  the  Na  concentration  in 
the  equilibrating  solution,  both  y    and  a    are  not  true  constants; 
hov;ever,  their  range  of  variation  can  be  determined o   For  the  experimental 
isotherm  with  Na"^  in  0.05  M  Ca(N02)2  over  the  0  -  200  ppm  Na   range, 
a  minimum  calculated  value  of  0.15  M  for  I  is  obtained  using  the  above 


47 


formula.   From  the  neaKiired  data  for  this  isotherm  a  naximum  value  of  I 
can  be  obtained .   Ifncn  the  highest  concentration  of  200  Mg/r.il  was 
used,  tlie  measured  equilibrium  values  were  S  =  22  Mg/g  and  C   =  178 
pg/ml .   Assuming  that  anions  are  excluded  completely  by  the  exchanger, 
the  concentration  of  ionic  species  in  the  solution  phase  was  Na 
0.0077  M,  Ca^"*"  =  Oo0505  M,  CI-  =  0.0057  M  and  NQ"  =  0.10  M  giving 


are  0.7A0  and  0.735,  respectively.   Corresponding  values  for  *a   are 

0.1^4  and  0.1'^2.   Thus,  the  values  of  y,,  and  ''a^  do  not  vary  appre- 
'  Na       Ca  J      fy 

ciably  between  the  maximum  and  minimum  I  values <.   Hov/ever,  it  is  clear 
that  with  decreasing  concentration  of  Ca(NO  )„  in  the  equilibrating 


/a        does  increase.   As  a  realistic  approximation,  an  average  value  of 


^a       =  0.1331.   Utilizing  these  values  and  the  experimental  K  value  of 

0.1324  for  this  isotherm,  a  value  of  K,  =  y,t  l^'^n      can  be  calculated. 

1    Na    Ca 

If  the  above  reasoning  and  assumptions  are  correct,  this  value  of  K, 
may  be  used  to  predict  the  experimental  K  values  for  the  adsorption 
isotherms  using  0.02  M  and  Oc005  M  Ca(NO  )„»   Similar  arguments  can  be 
advanced  for  the  equilibrium  isotherms  of  Li  „ 

V.'ith  the  above  approach,  K  values  for  the  isotherms  of  Na   in 
0.02  M  and  0.005  M  CaCNO^)^  were  predicted  using  a  calculated  K  value 
of  41.9  and  mean  I  values  of  0.064  M  and  0.020  M.   These  predicted  K 
values  were  0.201  and  0.347  which  compares  favorably  v;ith  the  measured 
values  of  0.2196  and  0.3770,  respectively.   Similar  computations  for 
Li"^  using  K,  =  70.18  and  mean  I  values  of  0.079  M  and  0.025  M, 


A8 


yielded  predicted  K  values  for  the  Li  adsorption  isotherms  in  0^02  M 
and  Oo05  M  CaiUO^)..  of  0„121  and  0.199  which  compared  reasonably  well 
with  the  neasured  values  of  Ocll2  and  0.153. 

Implicit  in  the  use  of  the  Debye-Huckel  theory  is  the  assumption 
that  the  activities  of  the  ions  in  the  exchanger  and  solution  phases 
are  predominantly  the  result  of  electrostatic  idn-ion  interactions. 
The  foregoing  theoretical  results  confirm  the  validity  of  this  assump- 
tion and  underscores  the  importance  of  such  interactions  in  ion- 
exchange  equilibriao 

45   2+ 
The  exchange  adsorption  isotherms  for   Ca   in  Oo075  M  and  0.05  M 

A5   2+ 
Ca(NO^)^  are  given  in  Figure  5„   These  were  linear  over  the   Ca   con- 
centrations used  and  illustrate  a  similar  decrease  in  the  K  value  with 
an  increase  in   Ca   concentration  in  the  equilibrating  solution. 

Assuming  that   Ca   and   Ca   are  indistinguishable,  then  Y'^^ca  ' 
y^^Q\-^   and  Y^"'ca  ^  'V^'^ca"   This  implies  that  there  are  no  differences  in 
their  physical  interactions  to  produce  selectivity „   With  the  previous 
theoretical  considerations  [equation  (33b)],  and  using  a  dagger  to  dis- 
tinguish between  properties  of   Ca   and   Ca 

^  =  i!^  =  1  (A9a) 

C-C   Y^ 


£-  =  ^  (49b) 

At  low  concentrations  of    Ca    both  C  and  C  can  be  regarded  as  con- 


~  '    C"*"  (A9c) 


S+  =  ^C"'  (49d) 
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E'iiplo}'iiig  the  same  approach  as  that  used  for  Na  and  Li  ,  it  is  possible 
to  compute  K  using  the  experinental  K  value  of  one  adsorption  isotherm. 
This  value  can  then  be  used  to  predict  the  K  value  of  other  isotherms. 

This  approach  provides  an  evaluation  of  the  validity  of  equation  (A9d). 

45   2+ 
The  experimental  K  value  of  the  isotherm  for   '  Ga    in  0.075  M  Ca(NO„) 

is  0.764  which  gives  K.  =  0.057.   The  predicted  K  value  of  the  isotherm 

for    Ca"   in  0.05  H  Ca(NO„)„  v/as  1.146  which  compares  favorably  with 

the  measured  value  of  1.041. 

45   2+ 
The  exchange  adsorption  isotherms  for   Ca   are  of  interest  in 

explaining  the  basis  for  the  effect  of  increasing   Ca   concentrations 

on  the  K  values.   Consider  specifically  the  exchanger  in  equilibrium 

with  0.05  M  Ca(NO„)„  solution.   If  a  Donnan  electrostatic  potential 

difference  is  set  up  across  the  interphase  region,  this  essentially 

equalizes  the  diffusion  rates  of  Ca    ions  into  and  out  of  the  exchanger 

2+ 
al lowing  the  continued  existence  of  a  higher  concentration  of  Ca    in 

the  exchanger  phase  to  maintain  electroneutrality.   Althougli  nacro- 

scopically  the  composition  of  either  phase  remains  fixed,  microscopic 

2+ 
excViange  of  Ca   ions  continues  to  occur  at  a  fixed  rate  across  the 

45   2+ 
phase  boundary.   If  some  fixed  amount  of    Ca    ions  is  introduced  in 

45   2+ 
the  solution,  it  is  reasonable  to  expect  that  at  equilibrium  the    Ca 

ions  entering  the  exchanger  depend  solely  on  their  relative  abundance 

to    Ca"   in  the  solution.   Increasing  the  solution  concentration  of 

Ca    to  0.075  M  would  reduce  the  relative  concentration  of   Ca    to 

Ca    on  the  exchanger  and  thus  the  macroscopic  adsorption  of    Ca 

is  decreased.   It  is  therefore  reasonable  to  assume  that  tlie  ratio  of 

""Ca"   in  the  exchanger  and    Ca    in  the  solution  is  inversely  propor- 
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Lionnl  to  the  concentration  of   Ca   in  the  solution.   Expressed 
qii-nititatively,  f^.05''^^  075  ~  0.075/0.05  which  is  wliat  was  concluded 
above  [equation  (49d) ]  using  a  different  approach.   The  foregoing  re- 
sults serve  to  give  credence  to  the  concept  of  ion  exchange  as  a 
Donnan-type  redistribution  of  ions  between  the  solution  and  exchanger 
phases. 

It  was  not  possible  to  determine  the  negative  adsorption  isotherms 
for  Cl~  in  batch  studies  because  the  increases  in  Cl~  concentration  were 
too  small  to  detect  above  random  variations  associated  with  the  specific- 
ion  electrode. 

Miscible  Displacement  Experiments  vjith  Na' 
Elution  curves  are  presented  in  Figures  6  through  14  for  a  series 
of  input  pulses  containing  Cl~,  HTO  and  Na"^  in  0.05  M,  0.02  M  or 
0.005  M  Ca(NO„)„, using  three  pore-water  velocities  ranging  from  1.5  to 
15  cm/hr  for  each  Ca(NO^)„  concentration.   These  displacements  were 

made  using  a  30.4-cm  long  column.   The  medium  was  packed  to  a  bulk 

3 

density  (p)  of  1.786  g/cm  and  had  a  saturated,  fractional  volumetric 

water  content,  9^,  of  0o34A.   Areas  under  the  breakthrough  curves  in 
Figures  6  through  14  were  determined  by  trapezoid  rule  integration  and 
in  each  case  indicated  complete  recovery  of  the  material  injectedo 
Thus,  complete  reversibility  of  the  mass  transfer  processes  was  achieved 
in  these  column  studieSo 

It  was  not  possible  to  maintain  identical  values  of  the  three 
pore-water  velocities  for  each  concentration  of  Ca(NO_)^.   Hovjever,  in 
no  case  did  the  measured  pore-water  velocity  deviate  by  more  than  ±  2%  of 
the  mean  values  (14o69,  7^29  and  1.52  cra/hr)  taken  over  the  three 
CaC^O^),,  concentrations..   Dispersion  coefficients  were  determined  from 


52 


:-  ./^ 


c   c 

•H 

C     'Ji 


O      G)      CO     r>     CD 


33 


r-i     . 

CM    O 


54 


II     .H 
01     II 


oa)cor>coio^^    (V    oi    ^ 


lO 

•H    -H 

c^ 

O     2 

H 

s:  -a 

u 

-ra   -u 

c   c 

ra  ^ 

-^     O 

D 

1       tH 

I— 1     CO 
O     O 

> 

+    i-l 

CO     O 


55 


CNlPl 


U     4J 

si3 


C 
>   -3- 


56 


nS 


lO 


If) 


SSt^ZIE*— 


O     X     z 

D         O         O 


J L 


^ 


J '  ' L 


oa)OOr>cDir)^(T)QJ^ 
o  ' 


u|u 


O  m 


Si-H 


o  en 


o    CO 


57 


r-)  o 
II     II 


O0)CpI>CDL07C0C>J 
~       *         *  *         *   O  * 


u|u 


53 


II  rH 

U 


C 

H   -a 


D  -a 

ex  -H 

0     rH 


o   a 

•H     S 

3  (y 


59 


00 


h- 


CM 


^^2a=s5=— x^ 


n    •  O 


b-nS^'^ 


J I I I I '       '       I 


o  * 


o 


o|u 


60 


61 


the  HTO  clution  curves  and  were  not  markedly  affected  by  the  Ca(NO_) 
concentration.   The  dispersion  coefficients  averaged  over  the  three 
Ca(NO^K  concentrations  and  corresponding  to  the  mean  pore-water 
velocity  values  given  above  were  0.453  +  0.025,  0.263  +  0,040,  and 

9 

0.0S2  +  0.007  cm^/hr,  respectively.   These  values  plotted  against  the 
average  velocity,  as  shown  in  Figure  15,  gave  a  linear  relationship. 
Extrapolation  to  zero  pore-water  velocity  gave  an  intercept  of  0.0458 
cm  /hr  (1.3  x  10  ^  era  /sec)  which  represents  the  diffusion  coefficient 
of  HTO.   Recent  values  reported  by  Mills  (51)  for  the  self-diffusion 

coefficient  of  HTO  are  1.724  ±  0.003  x  lO"^  cm^/sec  at  15  C  and  2.236  + 

-5    ^ 
0.004  X  10   cm^/sec  at  25  Cc   The  experiments  reported  in  this  study 

were  conducted  at  approximately  20  C.   A  linear  interpolation  yields  a 

-5   2 
value  of  1,88  x  10   cm  /sec  which  compares  favorably  with  the  experi- 
mental value  when  tortuosity  factors  are  considered. 

The  elution  curves  for  Cl~  appeared,  in  all  cases,  slightly  to  the 
left  of  those  for  HTO  confirming  the  expected  exclusion  of  Cl~  by  the 
exchanger.   Because  no  exclusion  isotherm  v/as  measured  for  Cl~,  the 
excluded  volumes  determined  by  the  magnitude  of  the  left-hand  shift 

of  the  Cl~  curve  from  C/C  =  0.5  at  V  =  V   for  these  nine  elution 
o  o 

curves  were  averaged.   The  excluded  volume  was  calculated  to  be  12.8 

3  _ 

cm  using  the  CI   data.   This  information  was  used  to  calculate  a 

"retardation"  coefficient  of  R  =  0.938  for  Cl~  from  the  relationship 

R  =  1  -  12.8/V  .   As  shown  in  the  figures  6  through  14  this  value  gave 

analytical  curves  which  described  the  elution  curves  for  Cl~. 

Analytical  curves  for  Na  ,  utilizing  R  values  (1  +  — )  calculated 

with  the  K  values  from  the  appropriate  isotherms,  described  the  elution 

curves  for  Na"*"  in  0.05  M  and  0.02  H  Ca(NO  )„  reasonably  v/ell.   These 
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results  indicate  that  equilibrium  was  instantaneous  and  followed  the 
experimental  isotherm,   The  somewliat  poorer  agreement  between  the 
analytical  and  experimental  curves  for  Na  in  0.02  M  Ca(NO„)-  was 
probably  the  result,  in  this  case  only,  of  not  using  the  sarae  batch 
of  0.02  M  Ca(NO^)„  solution  as  that  used  to  determine  the  adsorption 
isotherm.   Concentration  differences  between  the*  two  batches  would  ex- 
plain the  consistent  left  displacements  observed  in  this  case„   It  also 

underscores  the  sensitivity  of  the  system  to  changes  in  the  concen- 

2+ 
tration  of  Ca 

The  linear  equilibrium  model  failed  completely  to  describe  the 
elution  curves  for  Na   in  0.005  U   Ca(NO„)„o   This  was  unexpected  and  the 
displacements  in  0.05  M  and  0.005  M  Ca(NOT)„  were  repeated  using  a 
shorter  column  with  p  =  1«76A  g/cm  and  0^  =  0,344.   These  elution 
curves  are  presented  in  Figure  16  and  are  shoxm  to  behave  in  an  identi- 
cal manner  as  that  observed  for  the  longer  column.   This  rules  out  the 
possibility  that  the  observed  results  were  experimental  artifacts. 

Failure  of  the  equilibriuia  model  to  describe  experimental  data  is 
usually  construed  as  an  indication  of  non-equilibrium.   As  previously 


film  diffusion  would  constitute  the  main  resistance  to  mass  transfer. 
According  to  this  concept,  a  stagnant  fluid  film  is  thought  to  exist 
around  the  exchanger  particles  across  which  interphase  mass  transfer 
takes  place  by  molecular  diffusion.   The  thickness  of  this  film  would 
vary  inversely  with  fluid  velocity.   However,  no  marked  shifts  or 
changes  v/ere  observed  in  the  shapes  of  the  elution  curves  for  Na  or 
HTO  In  0.05  M,  0.02  H  or  0.005  M  CaCNO^)^  although  the  pore  velocity 
was  varied  by  almost  an  order  of  magnitude. 
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Whatever  thu  process  associated  with  the  observed  results  in 

Fi^jure  16,  it  clearly  depends  in  sone  fashion  upon  the  concentration  of 

2-1-  + 

Ca   in  the  system,  since  the  concentration  of  Na  was  held  constant  at 

80  -  85  ppm  in  all  cases.   The  theoretical  and  experimental  studies  of 

Helfferich  and  his  coworkers,  reviewed  previously,  indicated  that  the 

interdif fusion  coefficient  of  Na   in  either  phase  would  depend   upon 

+       2+ 
the  relative  concentrations  of  Na  and  Ca   in  that  phase  [equation 

+      2+ 
(26)].   Their  theory  showed  that  at  low  Na   to  Ca   ratios,  the  inter- 
diffusion  coefficient  is  close  to  the  diffusion  coefficient  of  Na  and 

2+ 
approaches  that  of  Ca   as  the  ratio  increases.   However,  it  was  in- 

+       2+ 
conceivable  that  the  diffusion  coefficients  of  Na  and  Ca   differed 

by  an  amount  large  enough  to  account  for  the  drastic  changes  in  the 

2+ 
elution  curves  observed  when  the  concentration  of  Ca   was  decreased 

from  0,02  N  to  0.005  M„   This  argument  coupled  with  the  observed  null 
effect  of  variations  in  pore  velocity  leads  to  the  conclusion  that  in- 
terdiffusion  in  films,  as  conceived  above,  was  not  the  process 
responsible  for  the  observed  behavior. 

The  possibility  existed  that  intra-particle  diffusion  may  be 
the  factor  responsible  for  the  results  presented  in  Figures  12  to  14 „ 
Hov/ever,  if  this  were  true,  it  was  difficult  to  explain  why  such 
results  appear  only  during  the  displacement  of  Ka  in  0.005  M  Ca(NO  ) 
and  not  with  the  two  higher  concentrations.   In  order  to  investigate 
the  resistance  to  intra-particle  mass  transfer  of  the  exchanger  par- 
ticles, unsulphonated,  spherical,  20  -  50  nesh,  macroporous  copolymer 
beads  were  used.   Their  mean  diameter  was  at  least  5  times  larger  than 
that  of  the  sulphonated  exchanger  particles  used  in  the  previous  ex- 
periments.  The  identical  tracer  solution  of  Na  in  0.005  II  CaCNO-)^ 
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was  displaced  through  a  20.5-crii  long  column  packed  with  these  headSo   The 
elution  curves  for  Cl~,  HTO  and  Na   are  shovm  dii  Figure  17.   These 
curves  showed  no  indication  of  diffusional  mass-transfer  processes  thus 
negating  the  possibility  that  intraparticle  diffusion  kinetics  were 
responsible  for  the  previously  observed  results^   The  results  in  Figure 
17  also  confirm  the  assumption  that  a  lov;  resist>ance  to  intraparticle 
mass  transfer  existed  for  the  raacroporous,  polystyrene  exchanger 
material s. 

From  the  foregoing  studies,  it  was  now  clear  that  diffusional 
mass  transfer  kinetics  was  not  the  major  factor  responsible  for  the 
failure  of  the  equilibrium  model  to  fit  the  elution  data  for  Na   in 
0.005  M  Ca(NO„)„.   Experiments  for  the  transport  of  Li  v/ere  next  used 
in  an  effort  to  gain  a  deeper  insight  and  provide  a  reasonable  explana- 
tion to  this  anomaly o 

Before  proceeding  onto  these  experiments,  displacements  of  two 
tracer  pulses  consisting  of  Na  ,  Cl~  and  HTO  in  0„05  M  and  0.02  M 
Ca(MO^)„  were  done  for  steady-state  unsaturated  flow  conditions  in  the 
same  column  used  in  the  foregoing  experiments »   The  column  was  kept 
under  a  constant  pressure  of  23  cm  of  water  using  the  method  described 
previously.   The  fractional  volumetric  water  content,  6  ,  was  reduced 
to  0.2740.   The  column  was  positioned  vertically  in  order  to  give  the 
highest  possible  flov^7  rate  without  incurring  large  water-content 
gradients  in  the  column„   The  flo\^7  rate  v;as  adjusted  so  that  a  constant 
pressure  head  of  5  -  6  cm  of  water  was  maintained  at  the  inlet „ 

The  elution  curves  for  these  displacements  are  presented  in 
Figures  18  and  19.   Areas  under  these  curves  by  the  trapezoid  rule  in- 
tegration indicated  tliat  the  material  injected  \;it1i  the  pulse  was 
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recovered  completely  in  the  effluent.   The  R  values  used  for  the 
saturated  flow  experiments  were  suitably  modified  to  take  into  account 
the  reduction  in  8  .   The  equilibrium  model  [equation  (45)]  v;as  used  to 
generate  the  analytical  curves  shown  in  Figures  18  and  19. 

Rather  poor  agreement  was  obtained  between  the  analytical  and 
experimental  curves  for  Na   in  0.05  M  Ca(NO„)„  and  was  worse  for  Na 
in  0,02  M  Ca(NO„)^.   The  general  tendency  was  a  displacement  of  the 
experimental  data  to  the  left  of  the  predicted  analytical  curve. 
It  is  clear  that  better  agreement  is  obtained  if  a  lov^er  R  value  were 
used.   It  is  possible  that  with  unsaturation  a  portion  of  the  porous 
matrix  becomes  either  Inaccessible  or  was  accessible  only  at  an  ex- 
tremely slow  rate.   Subtle  arguments  are  required  to  justify  a  shift  to 
the  left  of  the  predicted  equilibrium  curve  in  the  second  case.   These 
fortunately  were  discussed  elsewhere  (69,  70)  and  are  not  presented 
here  because  it  was  the  intention  of  this  study  to  make  only  cursor}'' 
observations  of  Na   transport  under  unsaturated  flow  conditions. 

M 1 s c ible  Displacement  Experiments  with  Li 
As  a  consequence  of  the  foregoing  results,  it  was  decided  to  ob- 
serve the  elution  behavior  of  pulse  inputs  of  Li   and  HTO  in  0.05  H, 
0.02  M,  and  0.005  M  Ca(NO  )   at  a  single  fluid  pore  velocity.   The 
experimental  and  analytical  curves  for  these  displacements  are  presented 
in  Figures  20,  21,  and  22,  respectively.   Areas  under  these  curves  by 
the  trapezoid  rule  indicated  complete  recovery  of  materials  injected 
with  the  pulse.   A  20.5-cm  long  column  packed  to  a  bulk  density  (p)  of 
1.764  g/cm  and  with  a  fractional  volumetric  water  content  (G  )  of  0.344, 
was  used  in  these  and  all  subsequent  displaceiicnts  involving  Li  „ 
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The  behavior  of  Li   in  0.05  M  Ca(NO  )„  agreed  with  the  theoretical 
curve,  but  the  linear  equilibrium  model  failed  to  describe  the  elution 
patterns  of  Li"*"  in  0.02  M  and  Oc005  M  CaCNO^)^.   A  mild  discrepancy  was 

apparent  in  the  front  portions  of  these  latter  curves  becoming  more 

2+ 
pronounced  as  the  Ca   concentration  was  decreasedo   In  order  to  obtain 

a  concentration  of  80  -  85  ppm  Li   in  the  tracer  solution,  it  was 

3 

necessary  to  add  12^5  cm  of  1  M  LiCl/liter  of  tracer  solution.   This 

resulted  in  density  differences  and  may  have  caused  the  observed  dis- 
crepancies which  were  similar  to  those  reported  by  Rose  and  Passioura 
(62). 

The  results  on  the  back  side  of  the  breakthrough  data  were  in 
disagreement  with  the  predicted  analytical  curves.   A  clue  to  the 
reason  for  this  anomalous  behavior  was  provided  by  the  observation  that 
the  deviations  commenced  at  approximately  one  pore  volume  after  change- 
over to  the  tracer  free  eluting  solution. 

As  discussed  above,  for  analytical  purposes,  it  was  necessary  to 
keep  the  concentration  of  Li   in  the  tracer  solutions  betv>'een  80  -  R5 
ppm.   This  concentration  of  Li  resulted  in  a  larger  numerica]  contri- 
bution to  the  ionic  strength  of  these  solutions  than  the  same  concen- 
tration of  Na  .   It  was,  therefore,  reasonable  to  assume  that 
differences  in  ionic  strengths  betvv^een  the  eluting  and  tracer  solutions 
were  the  reason  for  the  observed  behavior o   Figures  23  and  24  present 
elution  curves  for  the  displacement  of  the  same  tracer  solutions  of  Li 
in  0.02  M  and  0.005  M  Ca(N02)2,  but  with  the  ionic  strength  of  the  eluting 
solution  adjusted  to  match  that  of  the  Li   plus  CaCNO^)^  tracer  solu- 
tion.  The  observed  inflections  in  both  cases  wore  altered,  but  did  not 
disappear  completely;  however,  the  change  was  towards  better  agreement 
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with  the  calculated  curves.   In  order  to  pursue  this  idea  further  and 
to  confirm  that  the  concentration  of  Li  was  a  contributing  factor, 
displacements  were  conducted  using  50  -  55  ppm  Li   in  0.02  II  Ca(MO  )„ 
and  0.01  M  Ca(NO„)„.   These  elution  curves  are  shovvrn  in  Figure  25.   The 
analytical  curves  were  generated  using  a  value  of  D  obtained  from  a  pre- 
vious displacement  at  alnost  the  same  pore-water  velocity.   The  R  value 
for  the  displacement  of  Li  at  50  -  55  ppm  in  0.01  M  Ca(NO  )„  was  com- 
puted using  the  empirical  exponential  relationship  for  the  dependence 

2+ 
of  K  on  the  concentration  of  Ca    (Figure  4).   The  inflection  in  the 

elution  curve  for  the  displacement  of  Li  in  0.02  M  Ca(NO  )   was  less 
pronounced  than  for  the  displacement  of  Li  at  80  -  85  ppm.   The  in- 
flection in  the  elution  curve  for  the  displacement  in  0.01  M  Ca(KO„)„ 

was  more  pronounced  than  in  0.02  H  Ca(NO  )  ;  showing  similar  enhance- 

2+ 
iiieuL  with  decreasing  Ca   concentrations  as  those  observed  previously 

(Figures  21  and  22). 

The  foregoing  experiments  served  to  isolate  some  factors  associated 
with  the  observed  results,  but  did  not  provide  any  insight  into  the 
mechanism.   This  was  needed  to  explain  why  tlie  effect   appeared  only 
on  the  desorption  portions  of  the  elution  curves  and  why  the  inflections 
did  not  disappear  with  ionic  strength  corrections.   In  addition,  cal- 
culations utilizing  the  extended  Debye-HiJckel  formula  show  that  the 

24-  + 

activity  of  Ca   in  the  tracer  solutions  containing  80  -  85  ppm  Li 

in  0.05  M,  0.02  M  and  0.005  11  Ca(NO„)   were  respectively  2.4,  A. 8  and 
11.4%  less  than  the  corresponding  activities  in  the  pure  solutions. 

It  was  shovm  above  [equation  (48c)]  that  the  K  value  is  inversely  re- 

2+ 
lated  to  the  activity  of  Ca   ^   The  above  results  therefore  imply  a 

depression  in  adsorption  with  a  resulting  increase  instead  of  a  decrease 

in  the  Li   concentration  of  the  solution  phase. 
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These  facts  led  to  the  following  experiment  which  provided  a 
reasonable  explanation  for  both  the  observed  inflections  and  the 
anomalous  behavior  of  Na   in  0.005  H  CaCNO^)^"   The  solution  in  the 
column  was  displaced  with  a  tracer  solution  containing  80  -  85  ppm  of 
Li   in  0.005  M  Ca(NO„)„  until  the  concentration  of  Li   in  the  effluent 
was  equal  to  its  concentration  in  the  influent  s"olution.   This  im- 
plies that  equilibrium  was  achieved  at  all  points  within  the  system. 
This  solution  was  then  displaced  with  deionized  water  and  the  con- 
centrations of  all  components  in  the  effluent  were  monitored.   The 

2+ 
appcarance  of  Ca   and  Cl~  was  monitored  qualitatively  by  precipitation 


brown  ring  test.   At  exactly  one  pore  volume,  the  concentrations  of  all 
components  including  Li  fell  sharply  to  zero.   It  thus  became  clear 
tliat  the  electroneutrality  requirement  peculiar  to  ion  exchange  adsorp- 
tion processes  resulted  in  all  components  behaving  as  non-reactive 
solutes  once  equilibrium  was  achieved.   This  result  implies  that  an 

amount  of  Li   equivalent  to  -^  V  C  remained  adsorbed  by  the  exchanger. 
'  By  o  o 

The  adsorbed  Li  was  then  eluted  with  a  0.005  N  CaCNO^)^  solutiono 
The  elution  curve  showed  the  appearance  of  Li   shortly  before  one 
pore  volume,  rising  sharply  to  a  steady  maximum  concentration  and  then 
falling  sharply  to  zero  at  approximately  2.56  pore  volumeSo   This  be- 
havior is  analogous  to  that  observed  on  the  back  side  of  the  pulse 
elution  data  for  Li"^  in  0.005  U   CaCNO^)^  (Figure  22).   These  results 
are  presented  in  Figure  26,  and  indicate  that  the  exchange  adsorption 
for  Li   in  0.005  H  Ca(NO_)   was  distinct  from  its  exchange  desorption. 
If  the  same  K  value  characterized  both  processes,  the  elution  volume  for 
desorption  should  have  been  IL  .  =  1.79  pore  volumes  instead  of  observed 
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value  of  2.56.   This  indicates  that  a  higher  K  value  was  associated  v/ith 
the  desorption  process. 

The  samples  collected  from  the  desorption  experiment  were  tested 

2+  2+ 

for  Ca   with  Na-oxalate  solution.   The  tests  showed  that  Ca   appeared 

in  the  effluent  v.'ith  the  Li  .   It  was  also  observed  by  coraparing  the 

turbidity  developed  in  the  effluent  samples  with  a  standard  (0.005  H 

2+ 
Ca(NO^)^)  that  the  concentration  of  Ca   rose  sharply  to  a  maximum 

2+ 
concentration  which  was  less  than  0.005  11.   The  Ca   concentration 

remained  at  the  lower  level  until  the  elution  of  Li  was  complete, 

after  which  it  rose  to  0.005  M.   This  was  expected  because  electo- 

2+ 


utrality  necessitates  that  an  equivalent  amount  of  Ca   replace  the 


ne 

Li"^  desorbed.   These  results  indicate  that  the  desorption  of  Li' 

2+ 
occurs  with  a  concentration  of  Ca   in  the  solution  phase  less  than 

2+ 
0.005  M.   This  implies  a  lower  Ca   activity  and  therefore  a  higher  K 

value  [Equation  (48c) ]o 

The  concentration  of  Li   in  the  saturating  solution  was  0.0122  M  and 

from  the  exchange  desorption  curve  the  maximum  value  reached  was  0.0067 

K.      Applying  the  electroneutrality  requirement,  the  concentration  of 

+  2+ 

desorption  the  solution  phase  constitution  is  Li  =  0.0067  M,  Ca 

0.0016  M,  no"  =  0.01  M.   This  yields  a  value  for  the  ionic  strength  of 

C.0117  K.   With  this  value  the  extended  Debye-Hlickel  foriiiula  gives 

values  for  y, .  =  0.9012  and  >^^  =  0.0325.   The  empirical  formula  K  = 
'Li  Ca 

Y  ./70.1837  v'^   ,  used  previously  to  calculate  K  values  of  the  ad- 
sorption isotherms  for  Li  ,  yields  a  value  of  K  =  0.3952  which  gives 
—  =  2.03.   This  value  is  greater  than  the  experimental  value  of  1.56 
(Figure  26)  but  it  is  in  the  right  direction.   Observations  from  the 
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eluLion  curve  for  Li   in  0«02  M  Ca(NO„)„  given  in  Figure  21  suggest 
tliat  if  a  similar  experiracnt  as  that  described  above  were  conducted 
with  0.02  M  Ca(NO^)„,  a  maxiraum  concentration  of  Li  =  0.82  C  would 
result  for  the  desorption  process.   By  the  same  reasoning,  the  con- 
centration of  the  various  species  in  the  solution  :hase  would  be  Li   = 

2+ 
U.O]  M,  Ca   =  0.015  M,  NO  =  0.04  M  yielding  an  ionic  strength  of 

0.055  M.   Similar  calculations  yield  values  of  y   ^   =  0.8292,  »^   = 

0.0842,  a  predicted  K  value  of  0.1403,  and  ^  =  0.7205.   The  value  of 

■~  compares  favorably  with  the  observed  value  of  Oc.60  (Figure  21). 

The  inflections  observed  in  the  elution  curves  for  Li   (Figures 

21  through  25)  become  clear  in  light  of  the  above  results.   Wien  lower 

concentrations  of  Ca(NO„)„  were  used  in  the  tracer  pulse,  the  quantity 

of  Li  adsorbed  is  increased.   The  desorption  process  is  characterized 

by  a  K  value  that  is  dependent  on  the  activity  coefficient  of  Li  and 

2+ 
the  activity  of  Ca    in  the  desorbing  solution.   If  a  low  concentration 

of  Ca(NO„)„  was  used  in  this  solution,  the  equivalence  of  exchange  rc- 

2+ 
suits  in  a  reduced  value  for  the  activity  of  Ca   and  hence  a  corres- 
pondingly larger  K  valuCo   UTien  the  concentration  of  the  desorbing 
Ca(NO„)   solution  is  increased,  the  K  value  dc^creases.   In  addition,  at 

higher  concentrations  of  Ca(NO  )   ,  the  effect  of  the  adsorption  of 

2+ 
Ca   would  be  less  pronounced.   The  desorption  K  value  then  approaches 

more  closely  the  K  value  for  adsorption.   These  facts  explain  the  re- 
duced degree  of  inflection  when  the  concentration  of  the  eluting 
Ca(KO  )„  was  increased  to  match  the  ionic  strength  of  the  tracer  solu- 
tion, or  when  the  concentration  of  Li   in  the  latter  was  lowered.   In 
addition.  It  explains  why  the  inflections  did  not  disappear  after  the 
ionic  strength  adjusLment  and  why  they  were  not  observed  for  the 
displacement  of  Li   in  0.05  M  Ca(N0^)2. 
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'flif-^  foregoing  results  also  clarify  tlie  apparent  kinetic  effect 
ulvserved  for  displacements  of  Na   in  0.005  U   Ca(NO  )  .   Examination  of 
the  condition  under  which  the  data  in  Figures  11  through  14  v;ere  obtained 
reveals  that  because  of  large  R  values,  the  tracer  front  for  Na  ,  in 
all  cases,  had  not  reached  the  Cnvalue  in  the  effluent  at  the  time  of 
changeover  to  the  eluting  solution.   It  is  obvious  from  the  foregoing 
that  this  would  result  in  a  change  in  the  ionic  composition  of  the 
solution  phase  before  the  column  was  completely  equilibrated  v/ith  Na  . 
Thus,  instead  of  reaching  a  maximum  concentration  of  Co,  the  effluent 
would  reach  a  maximum  lower  than  Co.   This  idea  was  tested  by  a  step 
input  displacement  experiment  with  Na   in  0.005  M  Ca(NO„)  .   The  result 
given  in  Figure  27  shows  that  the  above  is  indeed  a  reasonable  ex- 
planation since  the  effluent  now  shows  none  of  the  previous  apparent 
kinetic  effect,   A  further  point  in  favor  of  the  correctness  of  this 
hypothesis  is  that  it  explains  why  the  apparent  kinetic  effect  was 
associated  with  only  the  front  portion  of  the  cuirves.   If  true 
kinetic  processes  were  operating  and  the  adsorption  isotherm  was  linear 
the  effect  would  have  been  symr^ietrical. 

In  light  of  the  foregoing,  the  data  for  Na"*"  in  0^02  H  Ca(NO  ) 
may  seem  anomalous  since  they  did  not  show  inflections.   It  is  likely 
that  inflections  such  as  those  for  Li   in  0.02  M  Ca(NO  )   may  have 
existed  but  were  overlooked  because  they  were  too  small. 

45  2+ 
Miscible  Displacement  Experiments  with    Ca 

45   2+ 

Two  displacements  of  pulse  inputs  were  conducted  with   Ca   in 

0.05  M  and  0.075  M  Ca(NO  )   to  verify  that  the  adsorption  isotherms  for 

^'5   2+ 

Ca   do  reflect  appropriately  in  their  transport  behavior^   Because 

of  large  R  values,  these  displacements  were  conducted  in  a  shorter 
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5.7-cm  long  column  packed  to  a  bulk  density  (p)  of  1.769  g/cm  and 
saturated  to  a  fractional  volumetric  water  content  (0  )  of  0.347. 

The  results  of  these  displacements  are  given  in  Figure  28.   Because 
of  the  shorter  column,  experimental  error  was  greater-   In  both  cases, 
if  higher  R  values  were  used,  closer  agreement  would  have  been  obtained 
between  the  experimental  and  analytical  curves. 

It  is  interesting  to  observe  tliat  for  the  displacement  in  0.05  M 
Ca(NO„)„  changeover  to  the  eluting  solution  before  the  tracer  front  had 
completely  appeared  in  the  effluent,  resulted  in  a  similar  disagreement 
between  the  analytical  and  experimental  curve  as  was  obtained  for  Na 
in  0,005  M  Ca(NO„)„.   Further,  the  use  of  a  much  larger  pulse  for  the 
displacement  in  0.075  M  Ca(NO„)„  produced  the  expected  inflection  at 
approximately  one  pore  volume  after  changeover  and  a  shift  to  the 
right  of  the  analytical  curve. 


+      2+ 
of  Na   to  Ca   on  the  in terdif fusion  coefficient  of  either  ion.   Exa- 
mination of  the  relevant  equation  of  Helfferlch  and  his  coworkers 
[equation  (2  6) ], reveals  that  in  the  case  of  isotypic  exchange,  the 

interdiffusion  coefficient  is  equal  to  the  diffusion  coefficient  of 

2+ 
Ca    and  independent  of  relative  concentrations.   Also,  increasing 

+      2+ 
the  relative  concentration  of  Na   to  Ca   should  result  in  the  inter- 

2+ 
diffusion  coefficient  approaching  the  diffusion  coefficient  of  Ca 

If  this  were  indeed  responsible  for  the  behavior  of  Na   in  0.005  M 

Ca(NO_)-j then  similar  effects  should  have  appeared  in  both  the  elutior 

curves  for   Ca   ,   This  however  was  not  observed. 
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CHAPTER  6 
SUMIIARY  AND  CONCLUSIONS 

Theoretical  approaches  to  describe  cation  exchange  processes  are 
fairly  recent  although  the  topic  spans  a  period  of  over  a  century. 
Despite  the  acknowledged  importance  of  this  process  in  soil  systems, 
there  are  only  a  few  references  which  deal  with  transport  of  ionic 
species.   In  this  study,  the  miscible  displacement  technique  was 
utilized  to  investigate  the  influence  of  exchange  adsorption  on  the 
transport  of  selected  inorganic  ions  in  a  porous  medium.   The  influence 
of  dispersion  and  exclusion  was  also  studied.   A  Ca-saturated  organic 
cation  exchanger  was  used  as  the  medium,  wliich  eliminated  the  possi- 
bility of  interactions  other  than  ion  exchange. 

Solution  and  exchanger  was  conceived  as  a  heterogeneous  system 
consisting  of  two  mixed  phases.   A  thermodynamic  treatment  based  on 
this  concept  predicted  the  exclusion  of  anions.   It  was  sho\im  that,  in 
general,  solution-phase  ions  common  to  those  initially  saturating  the 
exchanger  influence  the  adsorption  of  other  charged  species.   This  was 
adequately  verified  by  experimental  exchange  adsorption  isotherms  deter- 
mined for  Na"*"  and  Li"*"  in  0.05  11,  0.02  M  and  0.005  M  Ca(NO  )„  and  for 

_       _  _      2  ^ 

^5ca2-*-  in  0.075  M  and  0.05  U   Ca(NO  )  .   For  the  range  of  concentrations 
studied,  the  adsorption  isotherms  for  Li  and  Na"*"  were  linear  and 
their  slopes  increased  in  a  non-linear  fashion  with  decreasing  Ca'- 
concentration  in  the  equilibrium  solutions.   Similar  trends  wore  e'/i- 
dent  from  the  -'Ca   exchange  isotherms.   It  was  demonstrated  inductively 
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Mint  these  results  could  be  attributed  to  ion-ion  interactions  as 
conceived  and  quantified  by  the  Uebye-Huckel  theory* 

Miscible  displacement  experiments  involving  Na  ,  Li  ,    Ca" 
and  Cl~  were  conducted  to  examine  the  consequences  of  the  foregoing 
results,  and  to  determine  the  presence  of  kinetic  mass-transfer  pro- 
cesses.  All  experiments  were  performed  using  pulse  inputs  of  the 
tracer  solution.   Tritiated  water  (UTO)  was  present  in  all  tracer 
solutions  in  order  to  evaluate  dispersion.   An  asymptotic  solution 
to  the  cenvective-dispersion  transport  equation  for  a  reactive 
solute  (linear  adsorption  isotherm)  was  used  to  predict  the  experi- 
mental breakthrough  data. 

A  series  of  displacements  involving  pulses  of  Na"*",  HTO,  and  Cl~ 
in  0.05  M,  0.02  M  or  0.005  U   Ca(NO  )2  were  conducted  under  steady-stat€ 
water-saturated  flow  conditions  at  three  pore-water  velocities  ranging 
from  1.5  to  15  cm/hr.   Dispersion  coefficients  obtained  from  the  HTO 
elution  data  ;.'ere  not  influenced  by  Ca(NO^)   concentration  and  v;ere 
linearly  related  to  pore-water  velocity.   The  elution  curves  for  Cl~ 
were  all  displaced  to  the  left  of  those  for  HTO  confirming  its  exclu- 
sion by  the  exchanger.   Reasonable  agreement  was  obtained  betv.een  the 
computed  and  experimental  elution  curves  for  Na   in  0.05  II  and  0.02 
M  Ca(NO  )   at  all  pore-water  velocities  studied.   The  analytical 


solution,  however,  failed  to  describe  any  of  the  experiraental  elution 
curves  for  Na"^  in  0.005  H  Ca(NO  )2.   These  results  could  have  been 
used  to  support  the  concept  of  diffusional  mass-transfer  kinetics,  but 
it  \%fas  shown  experimentally  that  such  processes  were  not  involved. 

A  series  of  similar  steady-state,  water-saturated  column  experi- 
ments were  conducted  to  study  the  transjiort  bcliavior  of  Li   and  HTO 
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in  0.05  M,  0.02  M  and  0.005  H  Ca(NO„)„  at  pore-water  velocities  be- 
tween 7  and  8  cm/hr.   The  analytical  Kolution  described  the  front 
portions  of  the  experimental  elution  curves  reasonably  well  but, 
except  for  Li"^  in  0.05  U   Ca(NO^)„,  failed  to  describe  the  back  side 
of  the  elution  data.   Inflections  were  observed  on  the  back  side  of 
the  pulse  elution  datao   Further  experiments  indicated  that  the  sever- 
ity of  these  inflections  depended  upon  both  the  concentration  of 
Ca(NO.,)„  used  for  the  exchange  desorption  and  the  concentration  of 
Li  in  the  tracer  pulseo   These  observations  resulted  in  the  per- 
formance of  an  experiment  in  which  it  was  demonstrated  that  the  para- 
meter associated  v;ith  the  exchange  adsorption  of  Li   in  0.005  M 
Ca(NO„)„  was  different  from  its  counterpart  during  desorptiouo   It 
was  further  sho;m  by  inductive  reasoning  that  this  vzas  true  for  all 
CHL'es  and  v/as  a  direct  consequence  of  the  metathetical  nature  of  the 
sorption  process.   With  this  conslusion,  it  was  possible  to  clarify 
the  anomalous  and  unexplained  results  observed  in  the  displacement 
experiments  involving  Na  in  0.005  11  Ca(NO„)„o   These  concepts  v/ere 
further  supported  by  the  data  from  displacement  experiments  v/ith 
^^Ca""^  in  0  =  075  M  and  0.05  II  CaQM^)^. 

Displacement  experiments  were  conducted  to  obtain  a  preliminary 
description  of  the  transport  behavior  of  Na  ,  Cl~,  and  HTO  under  steady- 
state  unsaturated  water-flov/  conditions.   Txito  experiments  V7ere  con- 
ducted involving  pulse  inputs  containing  these  species  in  0.05  M  and 
0.02  11  Ca(NO^)„.   The  exclusion  of  Cl~  was  evident  in  both  experiments. 
The  analytical  solution  failed  to  describe  the  experimental  elution 


90 


curves  of  Ha'  in  both  cases.   These  cursory  observations  indicated 
that  unsaturation  resulted  in  portions  of  the  exchanger  phase  beconing 
either  inaccessible  or  accessible  only  by  diffusion  which  is  a  rate- 
controlled  process. 

The  foregoing  results,  viewed  in  their  entirety,  show  that  in 
general  the  exchange  adsorption  isotherm  for  a  cation  is  a  function 
of  the  activities  of  all  charged  species  in  both  the  solution  and 
exchanger  phases.   In  addition,  the  exchange  adsorption  process  is 
at  all  times  subject  to  the  electroneutrality  restraint,  and  this 
results  in  its  metathetical  character.   These  two  facts  are  funda- 
mental to  explaining  or  predicting  the  transport  behavior  of  cations 
in  reactive  media  including  soils.   The  physical  and  chemical  pro- 
perties of  soil  exchanger  materials  are  somewhat  different  from  those 
of  synthetic  exchangers.   The  applicability  of  the  concepts  and  re- 
lationships utilized  in  this  study,  to  ion  exchange  adsorption  and 
transport  in  soils  can  be  evaluated  only  by  further  investigation. 
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